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Preface 


A major  problem  confronting  some  modern  fighter-type  aircraft 
is  the  loss  of  control  due  to  roll-pitch  inertial  coupling. 

Some  of  the  aircraft  have  flight  control  systems  which,  when 
the  aircraft  is  performing  a roll,  sense  pitchup  and  therefore 
command  nose  down  stabilator  to  saturation.  When  this  occurs 
the  aircraft  loses  control.  The  Air  Force  Flight  Dynamics 
Laboratory  at  W.P.A.F.B.  is  interested  in  the  analysis  of 
a statically  unstable  fighter-type  aircraft  flying  at  high 
angles  of  attack.  Of  particular  interest  is  a comparison 
between  body  axis  rolls  and  velocity  axis  rolls  and  a possible 
method  of  analyzing  the  problem  of  uncontrollability  due  to 
the  flight  control  system. 

I wish  to  thank  my  thesis  advisor.  Dr.  John  D'Azzo  of 
the  Department  of  Electrical  Engineering  of  the  Air  Force 
Institute  of  Technology,  for  his  overall  guidance  throughout 
my  thesis.  I also  wish  to  thank  Captain  J.  Gary  Reid  of  the 
Department  of  Electrical  Engineering  of  the  Air  Force  Institute 
of  Technology  for  his  help  in  the  development  of  the  nonlinear 
equation  solver.  A special  thanks  is  due  to  Captain  James 
Silverthorn  of  the  Department  of  Aeronautics  and  Astronautics 
of  the  Air  Force  Institute  of  Technology  for  his  patience 
and  timely  efforts  in  helping  me  to  develop  a working  linearized 
model  of  the  aircraft. 
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Abstract 


A linearized  model  of  a fighter-type  aircraft  with  significant 
roll-pitch  inertial  coupling,  including  its  full  flight  control 
system, is  required  in  order  to  conduct  a comparative  analysis 
between  body  axis  rolls  and  velocity  (stability)  axis  rolls. 

The  stability  of  the  aircraft  is  checked  at  various  roll 
rates  for  both  body  axis  and  velocity  axis  rolls.  This  is 
done  by  examining  the  signs  of  the  eigenvalues  of  the  linearized 
model-positive  for  unstable  and  negative  for  stable.  It  is 
found  that  at  various  angles  of  attack  the  velocity  axis  rolls 
prove  to  be  at  least  as  stable  and,  in  most  cases,  more  stable 
than  body  axis  rolls.  The  stability  is  also  observable  for 
various  combinations  of  flight  control  systems. 

In  developing  a nonlinear  coupled  equation  solver,  a 
single  equation  with  known  solutions  is  considered  first. 

This  is  done  to  show  a simplified  version  of  what  the  nonlinear 
program  is  required  to  do.  Next,  a pair  of  nonlinear  coupled 
equations  is  analyzed.  The  development  of  the  program  for 
the  single  equation  case  proves  to  be  successful,  but  certain 
problems  arise  when  working  with  a pair  of  coupled  equations. 

This  thesis  provides  a good  foothold  on  a method  of  analysis 
known  as  Bifurcation  Analysis  and  Catastrophe  Theory  which 


can  be  used  to  solve  the  nonlinear  coupled  aircraft  equations. 

This  thesis  presents  some  of  the  problem  which  could  be  encountered. 


I also  would  like  to  thank  Mr.  Michael  Rise  and  Mr.  David 
Bowser  of  the  Control  Criteria  Branch  of  the  Air  Force  Flight 
Dynamics  Laboratory.  They  proposed  this  topic  and  were  very 
cooperative  when  I needed  special  information  about  the 
aircraft  model. 
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INVESTIGATION  OF  ROLL  PERFORMANCE  FOR  A HIGHLY  NONLINEAR 
STATICALLY  UNSTABLE  FIGHTER-TYPE  AIRCRAFT 

I . Int  roduct  ion 

Background 

As  the  design  of  fighter-type  aircraft  moves  toward  that 
of  smaller  and  thinner  wings,  the  nonlinear  aspects  play  an 
increasing  role  in  the  performance  of  the  aircraft  and  thus 
add  to  the  stability  problem.  The  equations  of  motion  of 
an  aircraft  (Ref  1:11)  a re 


l\.  - 

• • 

PI  - RJ  ♦ QR(I  - 1 ) - PQ  J . 

A Ai<  y <\  *. 

(i.n 

EM  - 

QI  ♦ riUVV  * (p2*R‘1Jx: 

• • 

(1.2) 

LN  - 

RI1  - PJx=  * WKV1,^  * Jxt 

(1.3) 

EFx 

- m(lVwq-VR) 

(1.4) 

EFy 

• 

- m(V*UR-WP) 

(1.5) 

FF, 

L 

• 

- m(W*VP-UQ) 

(1.0) 

It  can  be  seen  that  for  aircraft  with  larger  wings,  giving  a 
more  evenly  distributed  mass  over  the  wing  and  fuselage,  the 
moments  of  inertia  (I  ,If,  and  I.)  are  essentially  equal  to 
each  other  and  the  product  of  inertia  (.1  ) is  considered 

insignificant.  This  de-couples  the  six  equations  into  three 
longitudinal  and  three  lateral-directional  equations. 

In  some  fighters,  the  weight  distribution  is  primarily 
Ln  the  fuselage  since  the  wings  are  small  and  thin,  thus, 
the  differences  in  the  moments  of  Inertia  in  equations 
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1.1  thru  1.3  arc  sufficiently  large  and  they  cannot  he 
neglected.  These  differences  produce  what  is  known  as 
inertial  cross  - coupl i ng . If  a rolling  moment  is  introduced 
this  results  in  some  yawing  moment,  and  the  term  PR(I  -I  ) 
may  become  large  enough  to  cause  an  uncontrollable  pitching 
moment  (Ref  1:183). 

Another  important  coupling  term  is  the  product  of  inertia, 
J , which  appears  in  equations  1.1  thru  1.3.  For  some 

X z 

fighter- type  aircraft  the  product  of  inertia  is  large 
| « enough  to  play  a significant  role  in  the  analysis  of  the 

aircraft's  equations  of  motion.  This  adds  to  the  nonlinearity 
and  coupling  of  the  aircraft.  Therefore,  the  task  in  this 
thesis  is  to  analyze  the  effect  of  the  product  of  inertia 
and  all  moments  of  inertia  for  a modern  fighter-type  aircraft. 

; ( < 

Due  to  all  the  nonl incarit ics  m the  aircraft,  it  would 
also  be  of  major  interest  to  perform  a nonlinear  analysis  of 
the  aircraft.  The  Flight  Dynamics  Laboratory  has  developed 
a computer  program  (TIMEH)  which  gives  the  time  history  of 
an  aircraft  for  desired  inputs.  The  time  histories  have 
been  useful  in  showing  that  the  problem  of  instability  does 
exist  when  the  aircraft  is  performing  a roll  maneuver  at  a 
high  angle  of  attack;  however,  the  time  histories  do  not 
provide  any  information  as  to  where  the  problem  originates. 

A possible  approach  to  this  problem  is  known  as  Bifurcation 
Analysis  and  Catastrophe  Theory  Methodology  (BACTM)  (Ref  4). 
This  approach  considers  a set  of  nonlinear  coupled  equations 


j 


and  determines  their  solutions  as  a function  of  the  control 
parameters.  In  the  case  of  an  aircraft,  these  control 
parameters  would  be  the  pilot's  input  to  the  flaperon,  rudder, 
and  stabilator.  Since  the  actual  BACTM  program  is  not 
available  for  use,  the  first  step  is  to  develop  a comparable 
program,  using  the  basic  ideas  of  BACTM. 


Purpose 

The  purpose  of  this  thesis  is  to  analyze  the  problem  of 
roll  divergence  for  a statically  unstable  fighter-type  aircraft 
at  various  angles  of  attack.  Of  particular  interest  is  the 
performance  of  the  aircraft  at  high  angles  of  attack.  This 
divergence  is  caused  by  the  roll  pitch  coupling. 

The  problem  is  quite  large  and  requires  much  work.  This 
thesis  is  only  the  first  step  to  eventually  understanding  the 
nonlinear  phenomena  and  to  then  devise  corrective  control 
action.  This  first  step  contains  two  approaches: 

(1)  To  look  at  the  aircraft's  linearized  equations  of 
motion  and  determine  the  maximum  sustainable  roll  rate,  before 
instability  occurs,  for  a statically  unstable  fighter-type 
aircraft  at  various  angles  of  attack. 

In  this  approach  the  equilibrium  conditions  are  assumed, 
the  equations  are  linearized,  and  the  eigenvalues  are 
derived.  When  an  eigenvalue  crosses  the  imaginary  axis  into 
the  positive  real  side,  instability  occurs.  This  approach 
gives  some  very  good  information  but  is  rather  limited  and 
localized.  Since  the  equations  are  linearized  from  the 
start,  there  is  no  insight  to  the  nonlinearities 


global  analysis  of  the  aircraft. 

(2)  To  investigate  the  nonlinear  events  which  occur  in 
nonlinear  equations  by  using  the  methodology  of  bifurcation 
analysis  and  catastrophe  theory. 

This  approach  uses  the  nonlinear  equationCs)  to  solve  for  the 
equilibrium  conditions  and  linearizes  the  equation (s)  about 
those  conditions.  From  there,  the  eigenvalues  are  derived. 

This  is  a global  analysis  which  allows  the  whole  bifurcation 
surface  to  he  obtained.  It  also  gives  a deeper  understanding 
of  nonlinear  events  like  jump  phenomena  and  limit  cycle 
phenomena,  which  occurs  in  aircraft  at  high  angles  of 
attack. 

This  approach  is  a difficult  undertaking  so  this  thesis  is 
only  concerned  with  using  the  methodology  for  a simple  example. 


Organization 

The  first  part  of  this  thesis  presents  the  development 
and  use  of  the  linearized  aircraft  model.  An  introduction 
to  bifurcation  analysis  and  catastrophe  follows;  then  the 
development  of  a computer  program  to  incorporate  this  methodology. 

The  second  part  of  this  thesis  starts  with  the  analysis 
for  the  reduced  aircraft  model,  then  the  full  linearized 
model  including  full  flight  control  system.  The  nonlinear 
analysis  of  the  single  nonlinear  equation  follows.  Finally, 
an  attempt  at  a pair  of  coupled  nonlinear  equations  is  looked 
into. 


The  third  part  of  this  thesis  gives  the  results  of  the 
linear  analysis:  Hffect  of  aileron- rudder  interconnect, 
affect  of  angle  of  attack,  comparison  to  specifications, 
comparison  of  reduced  model  to  full  model,  and  maximum 
sustainable  roll  rate  for  both  body  axis  rolls  and  velocity 
axis  rolls  at  various  angles  of  attack.  This  is  followed 
hv  the  results  of  the  nonlinear  analysis. 

The  final  part  of  this  thesis  gives  the  conclusions  and 
recommendations  for  both  the  linear  analysis  and  the  nonlinear 
analysis . 

Linear  Model 

A set  of  stability  derivatives  and  a flight  control  system 
is  selected  from  a wind  tunnel  model  of  a statically  unstable 
fighter- type  aircraft.  The  Air  Force  Flight  Dynamics  Laboratory 
has  prepared  time  history  data  for  this  model.  At  high  angles 
of  attack  the  aircraft  becomes  uncontrollable,  thus  proving 
the  usefulness  of  the  stability  derivatives  and  PCS.  The 
control  surfaces  used  for  this  model  are  the  flapcron,  rudder, 
and  stabilator. 

Some  basic  assumptions  are  made  for  the  linearized  model 
of  the  aircraft  in  order  to  simplify  the  analysis:  Gravity  is 
ignored  which  eliminates  the  sine  and  cosine  dependence  of 
the  equations  of  motion.  Since  the  perturbations  are  small, 
the  second-order  terms  are  assumed  equal  to  zero.  The 
velocity  of  the  aircraft  is  assumed  to  be  constant  since  the 
primary  interest  is  in  the  steady-state  solutions. 
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Method  of  Linear Analysis 
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To  determine  the  stability  or  instability  of  the  aircraft,  it 
is  necessary  to  observe  the  eigenvalues  of  the  aircraft  model. 
This  can  be  done  by  developing  the  characteristic  equation 
from  the  linearized  model.  The  characteristic  equation  is 
a polynomial  equation  whose  factored  roots  are  the  eigenvalues 
of  the  model.  The  roll  rate  (P)  is  left  as  a variable  in  the 
characteristic  equation  so  that  the  migration  of  eigenvalues 
can  be  observed  for  increasing  values  of  roll  rate.  The 
flight  control  system  is  neglected  in  this  part  to  see  how 
the  aircraft  performs  without  it.  The  characteristic  equation 
is  limited  to  one  set  of  trim  conditions. 

To  determine  stability  or  instability  of  the  aircraft 
for  various  trim  conditions  it  is  necessary  to  observe  the 
eigenvalues  of  the  linearized  model  of  the  aircraft  with  its 
flight  control  system.  To  do  this,  the  linearized  model 
is  set  up  in  state  space  form 

• 

[A]X  - [B]X  (1.7) 

which  becomes 

X - [A* 1 B]X  (1.8) 

where  X is  an  nxl  state  vector  and  [A'*R]  is  an  nxn  plant 
coefficient  matrix.  The  matrix  A must  be  invertible.  The 
eigenvalues,  which  are  the  roots  of  the  plant  coefficient 
matrix,  are  obtained  from  [A  *B].  As  long  as  all  the  eigen- 
values have  negative  real  parts  the  aircraft  remains  stable. 
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but,  once  any  of  the  roots  become  positive,  the  mode  associated 
with  the  eigenvalue  becomes  unstable. 

Therefore,  it  is  important  to  be  able  to  distinguish  the 
eigenvalues  in  order  to  know  which  mode(s)  can  become  unstable. 
This  distinction  is  accomplished  by  obtaining  the  eigen- 
values for  the  model  without  feedback  and  also  by  obtaining  the 
eigenvectors  for  each  eigenvalue.  The  eigenvectors  are 
used  to  separate  the  longitudinal  eigenvalues  from  the 
lateral  eigenvalues  of  the  aircraft.  The  next  step  is  to 
obtain  the  eigenvalues  and  eigenvectors  of  the  model  with 
feedback.  Separating  the  longitudinal  eigenvalues  from  the 
lateral  eigenvalues  and  plotting  the  open-loop  roots  along 
with  the  closed- loop  roots  on  a graph,  the  modes  for  the  model 
with  feedback  are  then  identifiable.  A typical  plot  is 
shown  in  Appendix  D. 

To  obtain  the  maximum  sustained  roll  rate  for  a fighter- 
type  aircraft,  the  model  is  set  up  with  roll  rate  as  an 
input  variable.  Since  it  is  important  to  observe  both  body 
axis  rolls  and  velocity  axis  rolls,  the  states  P (roll  rate) 
and  R(yaw  rate)  are  considered  as  input  variables  to  the 
state  equations.  Other  elements  left  as  variable  parameters 
to  the  program  are  velocity,  altitude,  angle  of  attack,  and 
physical  parameters  of  the  aircraft  model.  With  this  set-up, 
the  eigenvalues  for  both  body  axis  and  velocity  axis  rolls 
are  obtained  at  various  angles  of  attack  and  various  roll  rates. 
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A comparison  of  stability  regions  for  body  axis  rolls  and 
velocity  axis  rolls  is  then  possible.  For  a given  angle  of 
attack  the  maximum  roll  rate  can  be  determined  hefore 
instability  occurred. 

The  model  of  the  aircraft  has  two  types  of  longitudinal 
flight  control  systems,  depending  on  whether  the  aircraft 
angle  of  attack  is  less  than  or  greater  than  20.4  degrees. 

It  is  of  interest  to  the  Air  Force  Flight  Dynamics  Laboratory 
to  check  the  roll  stability  of  the  aircraft  at  various  angles 
of  attack  using  each  of  the  two  flight  control  system 
configurations  separately.  This  is  easily  accomplished  since 
the  computer  program  for  the  linearized  analysis  (Appendix  E) 
is  designed  to  be  as  versatile  as  possible. 

Bifurcation  Analysis  and  Catastrophe  Theory 

This  section  gives  only  a brief  overview  of  bifurcation 
analysis  and  catastrophe  theory  which  is  essential  to  the 
understanding  of  the  nonlinear  analysis  in  this  thesis. 

A more  complete  explanation  of  these  theories  is  contained  in 
a report  by  Dr.  Mehra,  Mr.  Kessel,  and  Dr.  Carroll  titled 
"Global  Stability  and  Control  Anslvsis  of  Aircraft  at 
High  Angles  of  Attack",  (Ref  41. 

Consider  an  example  of  a single  nonlinear  equation. 


3^  - X3  ♦ (CjlX  ♦ C, 


(1.9) 


where  X is  the  state  and  C.  and  CT  are  the  control  parameters. 

A simplified  explanation  of  bifurcation  and  catastrophe  can  be 
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C2 


Fig.  1.1  X „ with  bifurcation 
eq 


in  terms  of  this  example  nonlinear  system.  Figure  1.1 

shows  the  solution  of  equation  1.9  for  a fixed  value  of 

and  varying  C2  from  -3  to  3.  The  points  on  the  curve  in 

Figure  1.1  marked  with  the  symbol  * are  points  of  bifurcation. 

This  is  where  the  slope  of  the  curve  is  infinite  and  occurs 

at  the  values  of  C2  * -2  and  2.  Between  these  two  points 

X has  3 solutions  -1  unstable  and  2 stable.  If  C,  = -3  is 
eq  Z 

the  starting  point  and  is  then  gradually  increased,  X 

64 

follows  along  branch  1 (a  stable  branch).  Once  C2  reaches 
a value  slightly  larger  than  2,  the  value  of  Xg^  jumps  rapidly 
to  a value  on  branch  2.  This  jumping  is  known  as  a catastrophe. 

Figure  1.2  shows  the  solution  to  equation  1.9  for  a 
different  fixed  value  of  while  varying  C-,  from  -3  to  3. 

Here  there  are  no  bifurcation  points  and  no  catastrophes. 

This  solution  has  a nice,  smooth  shape  and  is  a more  desirable 
solution  than  that  of  Figure  1.1. 


I 
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Fig.  1.2  Without  Bifurcation  Fig  1.3  Xgq  With  Inflection 

Point 


Figure  1.3  shows  the  solution  to  equation  1,9  for  • 0 
while  varying  C-,  from  -3  to  3.  Here  there  is  a point  where 
the  slope  is  infinite  but  the  control  parameter  does  not 
reverse  direction.  This  is  known  as  an  inflection  point. 

Figure  1.4  shows  what  is  known  as  a cusp  catastrophe  which 
is  representative  of  a third-order  equation.  For  other 
orders  of  equations,  there  are  other  shapes  and  other  types 
of  catastrophes  (Ref  4:24). 

Method  of  Nonlinear  Analysis 

Starting  with  a nonlinear  equation  or  set  of  simultaneous 
nonlinear  equations  and  a set  of  control  parameters,  a 
computer  program  is  needed  that  will  solve  all  of  the  equations, 
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Characteristic  Representation  of  X Showing  a Cusp 
Catastrophe . 


given  any  combination  of  control  parameters  whereby  the  control 
parameters  have  specified  upper  and  lower  limits.  A flow 
chart  which  represents  what  is  needed  in  the  program  is 
shown  in  figure  1.5  (Ref  4).  The  first  step  is  to  set  the 
control  parameters  of  the  equation(s)  and  all  the  other 
parameters  needed  in  the  program.  The  next  step  is  to  solve 


the  nonlinear  equation(s)  for  a given  pair  of  control 
parameters.  The  next  step  is  to  check  whether  the  equation(s) 
has  reached  a bifurcation  point,  then  increase  or  decrease 
the  value  of  the  control  variables  depending  upon  whether 
or  not  a bifurcation  point  has  been  reached.  Checking  for 
bifurcation  points  is  accomplished  by  itting  up  the 
Jacobian  matrix  of  the  system,  evaluated  at  the  equilibrium 
values,  and  obtaining  the  systems  eigenvalues.  A bifurcation 
point  occurs  when  an  eigenvalue  or  set  of  eigenvalues  passes 
through  the  imaginary  axis.  Continue  all  the  steps  above 
until  all  combinations  of  the  control  parameters  have  been 
considered.  Finally,  plot  the  solution  to  the  equations (s) 
vs  the  control  parameters. 

Validity  of  Analysis 

In  the  linearized  analysis  it  can  be  said  with  reasonable 
assurance  that  the  aircraft  equations,  stability  derivatives , 
and  flight  control  system  are  accurate  and  valid.  The  results 
show  stability  regions  for  both  body  axis  rolls  and  velocity 
axis  rolls  at  various  angles  of  attack,  and  are  very  informative. 
The  strongest  support  for  these  results  is  the  similar  trends 
found  by  the  Air  Force  Flight  Dynamics  Laboratory,  using 
a time  history  program  for  the  same  model  aircraft. 

For  the  case  of  a single  nonlinear  equation  analysis, 
the  program  that  shows  the  existence  of  bifurcation  points 
and  catastrophes  gives  results  which  have  an  accuracy  of  10 
The  accuracy  can  easily  be  increased  but  is  not  necessary  for 
the  demonstration  problem.  To  justify  the  solutions,  a 
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separate  analysis  ot*  the  nonlinear  equation  is  obtained  by 
using  a Texas  instruments  Sfl  hand  calculator.  This  method 
proves  to  be  time  consuming  but  it  allows  for  a comparison 
of  results.  To  further  strengthen  confidence  in  the  developed 
program,  the  general  shape  of  the  solution  to  the  equation 
is  similtar  to  that  of  an  equation  analysed  in  (Ref.  4). 
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1 1 . Linear  Analysis 
Reduced  Linearized  Model 

To  he  sure  that  the  aircraft  model  used  is  statically 
unstable,  it  is  necessary  to  analytically  check  the  stability 
of  the  model  without  a flight  control  system. 

The  equations  of  motion  of  an  aircraft  are 


EFx  - m (U+WQ-VR) 

• 

(2.1) 

EFv  - m(V+UR-WP) 

(2.21 

• 

FF,  - m (W+VP-UQ) 
z 

(2.51 

• • 

ZAl  - PIx-RJx:+QR( W‘P<3JX2 

(2.41 

• 

IAM  - QI  *PR(I  -UMr-RIJ 

y x 4 x z 

(2.51 

• • 

EAN  - RI.-PJX2+PQ(I  -IX)+QRJX. 

(2.61 

where  equations  2.1,  2.2.,  and  2.5  are  the  force  equations 
and  equations  2.4,  2.5,  and  2.6  are  the  moment  equations. 
Velocity  (U)  and  roll  rate  (P)  are  assumed  to  he  constant 
which  eliminates  equations  2.1  and  2.4.  It  is  assumed  that 
the  xy  plane  of  the  aircraft  is  the  plane  of  symmetry  and 
the  X-axis  is  selected  as  the  symmetric  axis.  This  allows 
the  product  of  inertia  (J  , 1 to  he  set  to  :ero.  Making 
the  following  assumptions 


where  aQ,  VQ,  Wq,  Qq , and  Rq  are  steady-state  values,  equations 
2.2  and  2.3  are  combined  with  their  respective  aerodynamic 
equations  (Ref  1).  Equations  2.4  and  2.5  are  then  substituted  into 
equations  2.2  and  2.3. 

The  Laplace  transform  of  the  equations  are 

(a1  s + aQ)  0(s)  + (b7  s4-  + b1  s ♦ bo)oi(s)  - 0 (2.8) 

(c,  s2  c1  s + cQ)  0 (s ) + (dj  s dQ)ot(s)  - 0 (2.9) 

where  the  a's,  b's,  c's,  and  d's  are  constants  which  are 
listed  in  Appendix  A.  The  characteristic  equation  produced 
from  this  model  is 

S4  * AjS5  * (A2P02  ♦ Aj)  S2  . (A4P02  * AS)S 

* (Vo4  * Vo2  * V ' 0 (2-10) 

where  A^,  A.,,...  Ag  are  constants  which  are  listed  in 
Appendix  A.  The  equation  is  purposely  set  up  as  a function 
of  steady-state  roll  rate  (P  ) in  order  to  study  the  effect 
that  roll  rate  has  on  the  aircraft.  The  full  development 
of  equation  2.10  is  in  Appendix  A. 

The  characteristic  equation  is  factored  to  obtain  four 
roots -two dutch  roll  roots  and  two  short  period  roots. 

These  roots  are  plotted  on  a root  locus  as  a function  of 
roll  rate.  This  plot  not  only  shows  the  static  stability 
of  the  aircraft  but  also  shows  the  effect  that  roll  rate  has 
on  the  aircraft  without  a flight  control  system  flying  at 


zero  angle  of  attack.  This  method  of  root  locus  analysis 
was  developed  by  Professor  John  Blakelock  (Ref  1).  This 
method  is  a helpful  tool  but  contains  too  many  assumptions 
for  this  particular  aircraft  model  and  is  limited  to  straight 
and  level  flights  only.  A more  versatile  approach  is  needed. 

Full  Linearized  Model 

The  linear  force  equations  and  angular  moment  equations  of 
motion  are  equations  2.1  thru  2.6.  Again  the  velocity 
(U)  is  assumed  constant.  Since  the  perturbations  are  small, 
the  second-order  terms  are  assumed  equal  to  zero. 

In  order  for  the  aircraft  model  to  be  analyzed  at 
various  angles  of  attack  and  various  roll  rates-both  in  body 
axis  and  in  velocity  axis  coordinates  - the  trimmed  (steady- 
state)  values  of  some  of  the  aircraft  states  are  left  as 
variables.  These  variables  include  velocity  (U),  angle  of 
attack  (a),  yaw  rate  (R) , and  roll  rate(P).  The  trimmed 
values  of  pitch  rate  CQ) , and  y-axis  velocity  (V),  are  set  to 
zero. 

Combining  equations  2.1  thru  2.6  with  their  aerodynamic 
force  and  moment  equations,  including  the  assumptions  made, 
the  following  linearized  aircraft  equations  are  obtained 

1 

(Ref  S): 

mUQB  + mUQr  - mU^ot  - mUoPQp  * qs  [Cy^B 

+ Cy  ( 7U7)p  + S/irH  r + Cy.  6f  + Cy5  5r]  (2,11) 

7p  O 7 r O 'Sf  5 


mU  a ♦ mU  P 0 -mU  q » qs[C  a + C • a 

O OO  O * * 2 Z 

a a 


+ Cz  (TD“)q  * cz,  M 
q o 5C 


(2.12) 


* ['’otIx'Iz1-2RoJx2]r  * tR0(Ix-Iz)*2P0J«]p 


qsc[cm  c * C (j§-)  q * C hqj-)a  * C Ss]  (2.1J) 
a q o a o p 


K-ijxz  * "Wy-Vx:”  O ■ S 

P 


* ci  (irf~  ) p + ci  (2Tr  ^ r * ce,  + ci  5f]  C2.14) 

p 0 r sr 


rIz  ‘ p Jxz  + ^o^y-V  + RoJxz]  q " qsb[Cn  6 

3 


* Cn„(2C}  P * CnJ  TlH  r + Cn,  **  * C 


P o 


f ^ V M (2*15) 


where  UQ,  PQ,  etc.  are  steady-state  values  and  p,  r,  etc.  are  small 
perturbations.  Appendix  B contains  a more  complete  derivation 
of  equations  2.10  thru  2.15.  The  values  of  the  stability  derivatives 
are  a function  of  angle  of  attack.  These  values  are  listed  in  the 
computer  program  EIGEN  in  Appendix  E and  are  given  only  for 
ot  ■ 0,  10,  20,  and  25  degrees.  For  the  purpose  of  this  thesis, 
these  values  of  angle  of  attack  are  sufficient. 
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Since  the  aircraft  model  is  statically  unstable,  it  is 
necessary  to  augment  the  aircraft  with  a flight  control  system. 


The  model  also  contains  many  nonlinear  terms  which  cannot 
be  ignored;  therefore,  both  a longitudinal  and  a lateral 
flight  control  system  are  included.  Figures  2.1  and  2.2  show 
a simplified  block  diagram  of  the  longitudinal  and  lateral 
flight  control  systems  used  in  this  thesis.  The  flight  control 
systems  were  supplied  by  the  Air  Force  Flight  Dynamics 
Laboratory. 

The  diagram  of  the  longitudinal  FCS  contains  two  different 
flight  control  systems,  the  use  of  which  depends  on  whether 
the  angle  of  attack  is  greater  than  or  less  than  20.4  degrees. 
The  lateral  FCS  contains  a component  known  as  an  aileron- 
rudder  interconnect  (ARI).  This  component  adds  to  the  non- 
linearity of  the  model. 

The  equations  for  these  flight  control  systems  are 
developed  in  terms  of  physical  variables  -XI  thru  X5  are 
the  variables  in  the  longitudinal  FCS  and  Z1  thru  24  are 
the  variables  in  the  lateral  FCS.  The  flaperon,  rudder,  and 
stabilator  displacements  are  also  included  as  physical 
variables.  The  development  and  listing  of  these  equations 
are  in  Appendix  C. 

The  next  step  is  to  combine  the  aircraft  equations  with 
the  longitudinal  and  lateral  FCS  and  to  arrange  them  in  the 
following  state  variable  form: 


b I 

i.  i 

L 


[AH*]  - [Aim  C2. 16) 


where  the  17  states  are 


The  elements  for  matricies  A and  B are  from  the  airplane 
equations  and  flight  control  system  and  are  listed  in  the 
computer  program  EIGEN  in  Appendix  E. 

The  computer  program  EIGEN  is  set  up  to  solve  for  the 
eigenvalues  of  the  following  system 

T-  [A*  1 B]  Y (2.32) 

which  are  obtained  from  equation  2.16. 

The  first  step  in  the  use  of  the  program  is  to  input  the 
desired  trim  conditions  of  the  aircraft.  The  inputs  include 
velocity,  altitude,  moments  of  inertia,  product  of  inertia, 
angle  of  attack,  roll  rate,  yaw  rate,  surface  area,  chord, 
wing  span,  air  pressure,  static  pressure  and  mass  of  the 
aircraft.  The  units  of  these  inputs  are  listed  in  Appendix  E. 

Once  this  is  done,  the  program  solves  for  the  eigenvalues 
and/or  eigenvectors  of  matrix  [A-  1 B ] . The  first  set  of  runs 
is  set  up  to  solve  for  the  eigenvalues  and  eigenvectors 
of  the  system,  both  with  and  without  a flight  control  system, 
for  a given  angle  of  attack  and  :ero  roll  rate.  This  permits 
the  determination  of  where  the  eigenvalues  originate. 

Without  feedback,  the  eigenvalues  of  the  flight  control 
system  are  easy  to  identify  since  the  FCS  is  set  up  in 
physical  variable  form.  This  leaves  five  eigenvalues 
unidentified.  At  this  point  the  eigenvectors  are  used  to 
identify  the  remaining  five  aircraft  eigenvalues:  Two 
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longitudinal  eigenvalues  (short  period  roots)  and  three 
lateral  eigenvalues  (two  dutch  roll  roots  and  one  rolling 
mode  root) . 

After  identifying  the  seventeen  eigenvalues  of  the  aircraft 
with  feedback  as  either  longitudinal  or  lateral,  the  next 
step  is  to  map  the  longitudinal  eigenvalues,  with  and  without 
feedback,  on  a graph.  The  same  is  done  with  the  lateral 
eigenvalues.  An  example  is  shown  in  Appendix  D.  This  mapping 
makes  it  easy  to  identify  each  of  the  eigenvalues  of  the 
aircraft  with  feedback.  A comparison  with  military  specifications 
is  made  for  the  short  period,  dutch  roll,  and  rolling  mode 
roots  in  order  to  prove  that  the  model  used  is  a valid  example 
of  a fighter-type  aircraft  (Ref  6). 

At  this  point,  a question  arises  as  to  whether  the  flight 
control  system  can  be  reduced,  thus  simplifying  the  analysis. 

Due  to  the  nonlinearities  in  the  aircraft  it  is  necessary 
to  retain  both  the  complete  longitudinal  and  lateral- 
directional  FCS.  One  possibility  is  the  removal  of  the 
aileron-rudder  interconnect  (ARI).  Without  the  ARI  one  of 
the  physical  states  (F4)  in  the  lateral  FCS  can  be  eliminated. 

An  analysis  of  the  short  period,  dutch  roll,  and  rolling 
mode  roots  of  the  aircraft  with  its  FCS  is  performed  both 
with  and  without  the  ARI.  This  is  done  at  0,  10,  20,  and  25 
degrees  angle  of  attack.  The  results  are  given  in  table  4.1 
with  the  conclusions  given  in  Chapter  V. 

The  next  step  is  to  study  the  effect  that  angle  of 


the  short  period,  dutch  roll,  an 
aircraft  with  IT'S  at  various  ang 
program  1: 1 GP.N  is  then  set  up  to 
of  the  aircraft,  with  PCS,  while 
Is  set  up  to  solve  for  the  eigen 
is  incremented  from  0 to  7 rad/s 
compare  the  full  model  with  the 
the  full  model  is  analysed  with 
attack,  zero  product  of  inertia, 
the  reduced  model  has  only  short 
these  are  t he  only  modes  which  v 
A stability  analysis  of  hod) 
axis  rolls  can  now  be  conducted, 
the  yaw  rate,  r,  is  set  to  zero 
incremented  from  0 to  7 rad/sec. 
to  roll  about  it's  own  X-axis, 
analyze  rolls  beyond  7 rad/sec  j 
Laboratory  is  not  interested  in 
point.  To  p e r f o r m v e l o c i l y ax i y 
incremented  from  0 to  7 rad/sec. 
set  as  a function  of  roll  rate 

r “ p t an  a 


d rolling  mode  roots  of  the 
les  of  attack.  The  computer 
solve  for  the  eigenvalues 
performing  rolls.  The  program 
values  while  the  roll  rate,  p, 
re.  The  first  step  is  to 
reduced  model.  For  comparison, 
inputs  of  zero  angle  of 
and  zero  feedback.  Since 
period  and  dutch  roll  roots, 
an  be  compared, 
axis  rolls  vs.  velocity 
To  perform  body  axis  rolls, 
while  the  roll  rate,  p,  is 
This  allows  the  aircraft 
It  is  not  necessary  to 
ince  the  A.  P.  Plight  Pvnamies 
an  analysis  beyond  that 
ro 11s,  the  ro 1 1 ra t e , p , is 
while  the  yaw  rate,  r,  is 
and  angle  of  attack: 

(2..VM 


The  combination  of  p and  r allows  the  aircraft  to  roll  about 
it's  velocity  vector. 
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The  program  EIGEN  is  then  used  to  scan  the  eigenvalues 
for  roll  rates  over  the  range  of  C to  7 rad/sec  about  the 
body  axis  and  the  velocity  axis,  and  to  record  the  areas  of 
stability  and  instability.  This  is  done  at  various  angles 
of  attack.  Of  primary  interest  is  the  case  of  high  angles  of  attack- 
20  and  25  degrees.  Plots  arc  made  with  body  axis  rolls, 
velocity  axis  rolls,  and  lines  of  instability  for  various 
angles  of  attack. 

The  final  point  of  interest  with  the  linear  analysis  is 
whether  it  is  necessary  to  have  two  longitudinal  flight  control 
systems,  depending  on  whether  the  angle  of  attack  is  less  than 
or  greater  than  20.4  degrees.  A plot  showing  lines  of 
instability  at  angles  of  attack  of  0,  10,  20,  and  25  degrees 
is  available  with  the  dual  longitudinal  flight  control 
systems.  A new  set  of  plots  (see  Chapter  IV)  is  then  run  at 
these  same  angles  of  attack,  but  only  one  of  the  longitudinal 
flight  control  systems  is  incorporated.  This  is  repeated 
for  the  other  flight  control  system. 


( 
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III.  Nonlinear  Analysis 


Single  Nonlinear  Equation 

In  order  to  use  the  ideas  of  bifurcation  analysis  3nd 
catastrophe  theory  a computer  program  is  needed  to  perform 
the  steps  shown  in  the  flow  chart  of  Figure  1.5.  To 
simplify  the  programming,  a single  nonlinear  equation  is 
analyzed  first.  The  nonlinear  equation  selected  (Ref  4:15) 
is 


dx  3 

- X ♦ CjX  + C2 


(3.1) 


where  X is  the  variable  and  and  are  control  parameters. 

The  first  step  is  to  set  the  control  parameters  to  some 

desired  initial  condition.  For  the  purpose  of  this  thesis, 

the  control  parameters  are  varied  from  -3  to  3;  therefore, 

the  initial  values  for  both  control  parameters  are  set  to  -3. 

The  next  step  is  to  solve  for  the  equilibrium  value  of 

X (X^) , given  the  values  of  both  control  parameters. 

An  AFIT  subroutine  NS01A  is  selected  to  solve  for  X . This 

eq 

subroutine  is  available  in  the  ASD  Computer  Center  library 
at  Wright-Patterson  AFB. 

In  order  to  check  for  a bifurcation  point  the  derivative 
of  equation  3.1  is  needed,  given  the  equilibrium  value  of  X 


3F 


3X“  + C, 

eq  1 


eq 
(3.2) 


eq 


where  g is  the  linearized  system  of  F.  If  g is  less  than 
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zero,  the  equilibrium  value  of  X is  stable.  If  g is  greater 


( than  zero,  the  equilibrium  value  of  X is  unstable.  If 

g equals  zero,  a bifurcation  point  is  encountered. 

The  digital  computer  program  is  set  up  with  two  loops- 

an  inner  loop  which,  increments  the  control  parameter  C2 

and  an  outer  loop  which  increments  the  control  parameter  C-^. 

In  the  inner  loop,  if  g is  not  equal  to  zero,  is  set  to 

C?+AC2,  where  aC2  is  a small  increment  of  C2.  If  g is 

equal  to  zero,  the  sign  of  AC  is  changed  and  then  C is 

« z 

set  to  C2+  AC2.  Everytime  a bifurcation  point  is  encountered, 
the  sign  of  AC2  is  changed.  For  example,  if  starts  at 
-3  and  is  increased  by  a small  amount  after  each  equilibrium 
value  of  X is  obtained,  C 2 is  increased  until  a bifurcation 
point  is  encountered.  At  this  point  C 2 would  start  decreasing 
until  another  bifurcation  point  is  met;  then  would  start 
increasing  again.  This  process  continues  until  C2  eventually 
reaches  it's  maximum  value. 

When  C2  completely  covers  the  range  of  -3  to  3,  the  outer 
loop  of  the  digital  program  is  engaged.  This  loop  resets 
C2  to  -3  and  increments  . Here  the  process  in  the  inner 
loop  starts  again.  Another  way  of  explaining  the  computer 
program  is  that  X is  obtained  as  a function  of  C2  with  C-^ 
set  to  a constant.  This  is  the  inner  loop.  All  the  outer 
loop  does  is  increase  the  value  of  and  returns  to  the 
inner  loop.  Once  the  range  of  -3  to  3 for  is  covered, 
the  process  of  solving  for  X is  completed. 

o 
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Fig.  3.1  Shifting  of  X at  Bifurcation  Point 


The  outputs  of  the  computer  program  include  a listing  of 

Xfi^  for  values  of  and  Cj,  the  residuals  of  X^,  and  the 

sum  of  the  squares  of  the  residuals  of  X . 

Since  the  subroutine  NS01A  solves  for  only  one  equilibrium 

value  of  X at  a time,  it  is  necessary  to  make  sure  that  the 

solution  of  equation  3.1  isn't  repeated  when  AC,  changes  sign. 

To  avoid  this,  the  next  guess  of  X is  changed  to  a value 

closer  to  the  next  nearest  value  of  X . An  illustration 

eq 

of  this  is  shown  in  Figure  (3.1)  where  * represents  a 

bifurcation  point  and  0 represents  a solution  of  equation 

3.1  already  known.  Once  C7  is  incremented,  the  subroutine 

NS01A  uses  the  present  value  of  Xg^  as  the  guess  to  compute 

the  next  value  of  X . At  C.*2  a bifurcation  point  is  reached 

eq  2 

and  therefore  C,  is  decremented.  In  order  to  prevent  X 
from  retracing  branch  (1),  it  is  necessary  to  adjust  the 


next  guess  of  X to  point  (■*•)  which  is  closer  to  branch  (2). 
e4 

This  allows  the  program  to  solve  for  X on  branch  (2). 

The  computer  program  which  performs  all  of  the  above  steps 
is  called  NONLIN  and  is  listed  in  Appendix  E. 

The  final  phase  of  the  program  is  to  plot  a series  of 
two-dimensional  arrays  showing  Xe^  vs  C.,  with  Cj  set  to  a 
desired  value.  These  plots  and  the  listing  of  X are 
compared  to  known  solutions  of  equation  3.1. 

Coupled  Nonlinear  Equations 

The  next  step  in  the  development  of  the  nonlinear  analysis 
technique  is  to  modify  NONLIN  to  solve  for  two  nonlinear 
coupled  equations.  The  equations  used  are 


dXl  3 

"3T  -V*  C1X2  * C> 


(3.3) 


dX2 

■ar 


X2  + C2X1  ♦ ci 


(3.4) 


The  part  of  NONLIN  which  checks  for  bifurcation  points  is  the 
only  part  which  needs  modification.  Once  the  equilibrium 
values  of  Xj  and  X-,  are  obtained  from  NS01A,  the  eigenvalues 
of  the  Jacobian  matrix  are  derived.  Since  the  inverse 
Jacobian  matrix  is  available  from  NS01A,  it  is  easy  to  solve 
for  the  eigenvalues  of  the  system.  Whenever  any  of  the 
eigenvalues  are  zero,  a bifurcation  point  is  reached  and 
the  sign  of  the  increment  of  the  inner  loop  control  parameter 
(C.,)  is  changed.  It  is  not  possible  to  know  whether  X.,  X,,  or 

X.  it 
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both  have  reached  a bifurcation  point.  This  problem  isn't 

considered  important  since  the  state  which  doesn't  hit  a 

bifurcation  point  will  only  backtrack  on  the  same  solution. 

Once  the  entire  range  of  C,  is  covered,  all  the  values  of 

and  X^^  have  been  derived.  The  two-dimensional 

plots  arc  set  up  to  give  the  solutions  of  vs  C-,  and 

X~  „ vs  for  a selected  value  of  C,  , of  C. , where  X, 

<-eq  2 11*  1 eq 

and  X1(i^  are  the  equilibrium  solutions  of  the  equations. 
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IV.  Results 


Linear  Analysis 

The  results  of  the  reduced  linearized  model  for  the  aircraft 
alone  are  shown  on  a graph  (Fig  4.1).  At  zero  roll  rate  the  short 
period  roots  are  real  (-0.57  and  +0.34).  And  the  dutch  roll 
roots  are  complex  conjugates  (-0.08  +_  j 1 . 0 6 ) . As  the  roll  rate  is 
increased  the  short  period  roots  move  away  from  the  imaginary 
axis  until  P = 0.6  rad/sec.  At  this  point  the  short  period 
roots  move  back  toward  the  imaginary  axis.  At  P = 1.23  rad/sec. 
the  aircraft  model  becomes  stable  because  both  short  period 
roots  are  now  negative.  At  P = 1.25  rad/sec.  the  short  period 
roots  change  from  two  real  roots  to  a pair  of  complex  conjugate 
roots.  As  P continues  increasing,  the  short  period  roots 
move  away  from  both  the  real  and  imaginary  axis.  The  dutch 
roll  roots  remain  complex  and  move  away  from  the  real  axis 
but  toward  the  imaginary  axis.  The  plot  of  Figure  4.1  contains 
the  roots  for  values  of  roll  rate  from  0 to  12  rad/sec. 

The  results  of  the  aircraft  with  the  FCS,  showing  the  effect 
of  variations  in  angle  of  attack  (AOA) , are  in  Table  (4-1).  This 
table  shows  the  results  for  the  model  both  with  and  without  the 
aileron  rudder  interconnect  (ARI).  The  short  period  roots  combine 
with  a pair  of  longitudinal  FCS  roots  to  form  two  pairs  of 
complex  conjugate  roots.  Their  damping  ratios  and  natural 
frequencies  remain  about  the  same  for  0,  10,  and  20  degrees 
AOA.  At  25  degrees  AOA,  the  short  period  roots  break  away 
from  the  longitudinal  FCS  and  form  a single  pair  of  complex 
conjugate  roots  and  two  real  roots.  The  rolling  mode  root 
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(REAL  AXIS) 


combines  with  a longitudinal  PCS  root  to  form  a pair  of 
complex  conjugate  roots.  At  10,  20,  and  25  degrees  AOA  the 
two  complex  roots  break  into  one  real  PCS  root  and  one 
real  rolling  mode  root. 

The  results  for  the  aircraft  model  with  the  ARI  are 
basically  the  same  as  for  the  aircraft  without  the  ARI  as 
shown  in  Table  (4.1).  The  differences  show  up  in  the  roots 
of  the  lateral  FCS.  Without  the  ARI,  one  of  the  eigenvalues 
from  the  lateral  FCS  is  unstable  for  a straight  and  level 
flight.  With  the  ARI  connected,  the  lateral  FCS  becomes 
stable.  From  this  point  on,  the  results  presented  in  this 
report  include  the  aircraft  model  with  the  aileron- rudder 
interconnect . 

For  a category  A flight  phase  with  level  one  flying 

quality,  the  minimum  short  period  dampling  ratio  Cp  _) 

sp 

is  0.35,  the  minimum  dutch  roll  damping  ratio  (p  is 
0.19,  and  the  minimum  dutch  roll  natural  frequency  (co^  ) 
is  1.0  rad/sec.  (Ref  6).  These  are  the  specifications  for 
a fighter  type  aircraft  under  highly  maneuverable  flight 
conditions.  For  all  angles  of  attack  the  short  period  damping 
ratios  meet  the  specifications.  The  dutch  roll  damping  ratios 
also  meet  the  specifications.  The  dutch  roll  natural  frequencies 
fall  short  of  the  specifications  at  0 and  25  degrees  AOA 
but  are  within  specifications  at  10  and  25  degrees  AOA. 
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The  results  of  the  reduced  linearized  model  are  then 
compared  with  the  results  of  the  full  linearized  model. 

It  is  found  that  both  models  produce  a pair  of  complex 
conjugate  dutch  roll  roots  (.*0.084  +_  j 1.06  for  the  reduced 
model  and  -1.32  +_  j 1.06  for  the  full  model)  and  a pair 
of  real  short  period  roots  (-0.567  and  0.334  for  the  reduced 
model  and  -0.664  and  0.176  for  the  full  model). 

To  show  the  results  for  body  axis  rolls  vs  velocity 
axis  rolls  the  regions  of  instability  are  needed.  Figure 
4.2  shows  how  these  lines  relate  to  each  other  for  various 
angles  of  attack.  The  area  to  the  right  of  the  lines 
is  where  the  rolling  mode  root  is  unstable.  This  is  for  the 
range  of  0 to  7 rad/sec.  roll  rate. 

Figure  4.3  shows  the  results  of  body  axis  rolls  vs  velocity 
axis  rolls  at  0,  10,  20,  and  25  degrees  AOA.  This  is  with 
the  dual  longitudinal  FCS.  At  0 degrees  AOA,  a body  axis 
roll  is  the  same  as  a velocity  axis  roll.  This  is  apparent 
from  equation  3.1.  Figure  4.4  shows  the  same  results  using 
the  longitudinal  FCS  for  angles  of  attack  less  than  20.4  degrees. 
Figure  4.5  shows  the  same  results  using  the  longitudinal 
FCS  for  angles  of  attack  greater  than  20.4  degrees. 

The  results  contained  in  Figures  4.3  and  4.4  are  compiled 
in  Table  4.2.  At  0 degrees  AOA  the  aircraft  is  unstable  at 
1.8  rad/sec.  roll  rate  for  the  dual  FCS.  At  10  degrees  AOA 
the  aircraft  is  unstable  at  2.0  rad/sec  roll  rate  for  a body 
axis  roll  and  at  2.14  rad/sec  roll  rate  for  a velocity  axis 
roll.  At  20  degrees  AOA,  the  aircraft  is  unstable  at  2.4 
rad/sec  roll  rate  for  a body  axis  roll  and  at  3.29  rad/sec 
roll  rate  for  a velocity  axis  roll.  At  25  degrees  AOA,  the 


R( RAD/SEC ) 

.00  cP.OO  4.00  8.00 

a w J.  i a 


Fig.  4.3  Body  Axis  vs  Velocity  Axis  Rolls  with  Dual  Long.  FCS. 
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Fig.  4.5  Body  Axis  vs  Velocity  Axis  Rolls  with  High  AOA  Long.  FCS 


ROLL  RATE  WHERE  AIRCRAFT  GOES  UNSTABLE  (RAD/SEC.) 
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1.8 
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10 
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2.14 

1.8 

1.86 

2.0 

2.14 

20 

2.4 
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1.8 

2.14 

2.4 

3.29 

25 

2.2 

>7.0 

2.2 

>7.0 

2.2 

> 7.0 

Table  4.2  Occurance  of  Instability  for  Aircraft  with  Various 
FCS's. 


aircraft  is  unstable  at  2.2  rad/sec  roll  rate  for  a body 
axis  roll  and  remains  stable  to  7.9  rad/sec  roll  rate  for 
a velocity  axis  roll.  In  all  cases  the  aircraft  remains 
stable  for  larger  velocity  axis  rolls  than  for  body  axis 
rolls.  As  the  AOA  increases,  the  aircraft  remains  stable 
for  larger  values  of  roll  rate  for  velocity  axis  rolls. 

At  25  degrees  AOA  the  aircraft  is  stable  up  to  8.6  rad/sec. 
roll  rate  for  a velocity  axis  roll.  These  same  trends  were 
found  by  the  Air  Force  Flight  Dynamics  Laboratory  in  their 
time  history  evaluations  of  the  same  model. 

A comparison  of  the  three  flight  control  system  configurations 
show  that  the  high  AOA  longitudinal  FCS  gives  poorer  results 
for  body  axis  and  velocity  axis  rolls  at  0,  10,  and  20  degrees 
AOA  when  compared  with  the  dual  longitudinal  FCS.  The  results 
also  show  that  the  low  AOA  longitudinal  FCS  produces  instabilities 
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at  smaller  roll  rates  than  the  dual  longitudinal  FCS  when 
at  25  degrees  AOA.  These  results  are  shown  in  Table  4.2. 


Nonlinear  Analysis 

The  results  of  the  single  nonlinear  equation  are  shown  in 
Figure  4.6  which  is  produced  by  use  of  the  computer 
program  NONLIN.  Bifurcation  points  are  encountered  for 
Cj  < 0 and  cusp  catastrophes  occur.  These  results  are 
shown  to  be  correct  when  compared  with  the  results  obtained  with 
a Texas  Instruments  59  hand  calculator.  The  characteristic 
shape  of  the  curve  of  X vs  C is  similiar  to  that  shown  in  a 
report  by  Dr.  Mehra,  Mr.  Kessel,  and  Dr.  Carroll  (Ref  4:17). 

When  the  two  coupled  equations  are  tried,  a bifurcation 
point  is  never  encountered.  While  incrementing  the  control 
parameter  C7 , an  eigenvalue  does  cross  the  imaginary  axis 
which  is  an  indication  that  X^,  or  both  X^  and  X., 
has  reached  a bifurcation  point.  The  data  shows  that  there  is  no 
indication  of  a bifurcation  point.  There  also  is  no  indication 
of  an  inflection  point. 
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V.  Conclusions  and  Recommendations 


Linear  Analysis  Conclusions  I 

Since  the  lateral  FCS  is  statically  unstable  without  the  j 

aileron  rudder  interconnect  CARI)  it  is  obvious  that  the 
ARI  is  a necessary  part  of  the  model’s  FCS  and  cannot  be 

excluded.  U 

A comparison  of  the  model's  roots  with  specifications 
shows  that  the  aircraft  model  used  in  this  analysis  is 
representative  of  a fighter-type  aircraft. 

When  the  reduced  model  is  compared  to  the  full  model, 
it  is  observed  that  the  short  period  roots  and  dutch  roll 
roots  are  closely  related  but  different.  The  dutch  roll 

roots  are  different,  primarily  because  the  aerodynamics  of  the  U 

M-moment  equation  and  the  Z-force  equation  are  neglected  in 
the  reduced  model.  The  short  period  roots  are  different, 
primarily  because  the  aerodynamics  of  the  N-moment  equation 

l 

are  neglected. 

In  the  comparison  of  body  axis  rolls  vs  velocity  axis 
rolls,  the  aircraft  remains  stable  for  larger  roll  rates 
when  it  rolls  about  its  velocity  vector.  The  comparison  of 
the  three  sets  of  longitudinal  flight  control  systems  show 
the  necessity  for  the  dual  longitudinal  FCS.  This  is  not 
to  say  that  there  doesn't  exist  a single  FCS  which  gives 
better  roll  rate  performance.  The  results  show  that  the 
combination  of  the  two  longitudinal  flight  control  systems 

performs  better  over  the  range  of  0 to  25  degrees  AOA  than 
when  just  one  of  the  longitudinal  flight  control  systems  is  used. 
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Linear  Analvsis  Recommendations 


The  linear  analysis  in  this  thesis  has  produced  an 
acceptable  linear  model  of  a statically  unstable  fighter- 
type  aircraft.  This  thesis  has  also  shown  that  the  aircraft 
becomes  unstable  when  performing  rolls.  The  elimination 
of  the  problem  has  not  been  looked  into. 

One  possible  procedure  would  be  to  set  up  a computer 
program  that  takes  the  state  space  model  of  the  aircraft 
and  produces  time  responses  for  any  state  to  a desired 
input.  This  has  been  done  in  a computer  program  called  TOTAL 
which  is  available  in  the  Air  Force  Institute  of  Technology 
program  library.  The  dimension  statements  in  TOTAL  presently 
limit  the  model  to  no  more  than  ten  states.  Since  this  model 
has  seventeen  states,  a change  of  the  dimension  statements 
programming  is  needed. 

Once  this  is  done,  a study  of  the  influence  that  various 
feedback  components  have  on  the  states  of  the  model  can  be 
performed.  This  would  possibly  help  find  the  cause  of  the 
loss  of  control.  The  final  step  would  be  to  eliminate  the 
loss  of  control  without  depreciating  the  performance  of  the 
aircraft . 

Nonlinear  Analysis  Conclusions 

The  analysis  of  the  single  nonlinear  equation  has  provided 
a good  initial  step  in  the  development  of  a computer  program 
for  applying  bifurcation  analysis  and  catastrophe  theory  to 
a control  problem.  The  coupled  equations  used  proved  to  be 
unsatisfactory  for  this  analysis.  A more  indepth  study 
is  required  at  this  point. 
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Nonlinear  Analysis  Recommendations 


This  study  has  laid  the  ground  work  for  using  bifuraction 
analysis  and  catastrophe  theory  to  analyze  nonlinear  equations. 

It  is  possible  that  this  same  approach  could  be  used  to  analyze 
the  nonlinear  aircraft  equations  of  motion.  What  is  needed 
is  a computer  program  that  uses  the  same  principles  outlined 
in  Figure  1.5  and  is  capable  of  solving  the  aircraft 
equations.  From  here,  a study  of  the  influence  that  various 
feedback  components  have  on  the  solutions  to  the  aircraft 
equations  could  be  performed.  Finally,  determine  where  the  loss 
of  control  is  coming  from  while  performing  rolls  and  then 
eliminate  the  problem  without  depreciating  the'  overall  performance 
of  the  aircraft. 
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APPENDIX  A 

Development  of  Reduced  Linearized  Equation 


The  equations  of  linear  motion  and  the  equations  of 
angular  motion  are: 


IAFx  = m(U+WQ-VR) 

(A-l) 

ZAFy  = m (V+UR-WP) 

(A-2) 

ZAFz  = m (W+VP-UQ) 

(A-3) 

ZAL  " PIx  ' k Jxz  + 

qr(iz  - y 

- PQ  Jxz 

(A-4) 

ZAM  = $1  + PR(IX  - 

V * CP2  - 

R2)J 

1 xz 

(A-5) 

ZAN  - RIZ  - P Jxz  - 

PQdy  - lx) 

+ QR  J 
x xz 

(A-6) 

Assuming  constant  velocity  (U)  eliminates  the  X-force 


equation  (A-l).  Assuming  constant  roll  rate  (P  ) eliminates 


the  AL  moment  equation  (A-4).  This  assumption  was  made 
because  the  transient  value  of  roll  rate  was  not-needed 


for  the  stability  analysis.  Assuming  VQ  = Wq  = Qq  = Rq  = 0 , 


then  Equations  A-2,  A-3,  A-5,  and  A-6  became 


ZAF  = m(v  + Ur  - wP.) 
y o 0 


ZAF  = m(w  + vP  - U q) 

Z O 0 


EAM  = qly  * P0r(Ix  - Iz)  * (Po2  - r2)  Jxz 


IAN  - r I,  ♦ Poq(Iy  - Ix)  ♦ qr  Jxz 


(A-  7) 
(A-8) 
(A-9) 
(A- 10) 


where  Uq  and  PQ  are  equilibrium  values  and  v,  r,  w,  and  q are 


small  perturbations  from  the  equilibrium  values 
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Assuming  that  the  aircraft  has  the  XY  plane  as  the  plane  of 

symmetry  and  selecting  the  X-axis  as  the  symmetry  axis,  then 

J = 0.  Assuming  EAM  = EAN  = EAF  = EAF  = 0,  equations 
XZ  X z 

A- 7 thru  A- 10  became 


EAM  = qly  + Por(Ix  - Iz) 

FAN  - flz  ♦ Poqdy  - Ix) 

EAF  = § + r-  Pa  = 0 
y o 

EAF  = a + P(3-q  = 0 
2 0 1 


where  S = v/U  and  a ~ fi^a  ~ w/U  . 

o o 


(A- 11) 


(A- 12) 


(A- 13) 
CA  -14) 


The  aerodynamic  moments,  including  all  the  stability 
derivatives  for  the  m and  n moment  equations,  are  (Ref  1:21). 


-C  U - 
mu 


xz  V 
<t> 


C a - C a - 757^—  C 0 
ma  2Uo  mq 


(A* 15) 


k c 
2U  n 
o p 


2ir  C ♦ - Cn/ 


(A- 16) 


Equating  equation  A- 11  and  A- 15  and  letting  q = V yields: 


— ^ — q ” TTi — C q - — C ct  - C ot  + 

Sqc  ^ "q  2Uo  n« 


I - I 

— 2.  P r = 0 (A- 17) 

Sqc  ° 


Equating  equation  A-12  and  A-16  and  letting  r = V and  Pq  = 
yields : 


b Cl  - I ) 

■£-  C r - C 0 + — ^ 5- 

2Uo  nr  ' nB  ' Sqb 


rn-  c p 

2Uo  np  ° 


(A- 18) 


Solving  equations  A- 13  and  A- 14  for  r and  q yields: 


r = P a-5 
o 


q = P B+a 
o 


(A- 1 9) 


(A- 20) 


and  differentiating,  results  in 


r = P ot-B 
o 


q = P S+a 
o 


(A-21) 


(A- 22) 


Substituting  equations  A-19  thru  A-22  into  equations  A-17 
and  A-18  yields: 


— (P  B+a)  - C (P  B+a)  - - C a - C a 

Sqc  0 2Uo  mq  0 “Uo  ma  ma 


(I  -I  ) 

Sqc"^  Po<V‘S>  " 0 


(A- 23) 


c„  cp0.-S)  - c^s  . P0(P0s*i) 


th~  C P 
2Uo  np  0 


(A- 24) 


Regrouping  and  taking  the  Laplace  transform  of  equations  A-23 
and  A-24  vields: 


(I  -I  + I )P 

v y x ZJ  0 

Sqc 


s - c p 

2U  m o 
o q 


3(s) 


-4-  s2  - (Cm  ♦ C ) s + (Ix-  '_Iz)P-?  - C 

p zu  m m p 

Sqc  o q a Sqc 


m 


a(s)  = 0 


(A-25) 


1z  2 b r e 

^ ^ V ' 


(I  + I 
z y 


(I  - I ) P ‘ 
y x o 

Sqb 


+ C 


n. 


6 (s) 


Sqb 


UP0  b 

x 0 S + C Prt 

2U0  nr  0 


a (s)  = 0 (A- 26) 


Equations  A-25  and  A- 26  have  the  form 


(axs  + ao)8(s)  + (b2sZ  + b:s  + b0)a(s) 


(c7sz  + c.s  + c ) 8 (s)  + (d,  s + dja(s)  = 0 


(A- 27) 


(A- 28) 


and  tha  characteristic  equation  is  fourth-order  and  is  a function 

of  P . The  stability  derivatives  and  aircraft  dynamic  terms 
o 

are  listed  in  the  computer  program  EIGEN  (Appendix  E) . When 
these  values  are  substituted  into  equations  A-25  and  A-26, 
the  characteristic  equation  becomes 

s4  + .382498s3  + (1.6961P  2 + .11241)s2 

o 


+ ( . 403345Pq  + . 017979) s 

+ ( . 718984Pq4  - . 056358Pq2  + .010219)  = 0 


(A- 29) 
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APPENDIX  B 

Development  of  Full  Aircraft  Linear  Equations 


The  equations  of  linear  motion  and  the  equations  of 
angular  motion  are: 


EAF  * 

X 

m (U+WQ- VR) 

(B-l) 

EAF  ■ 

y 

m(^'+UR-WP) 

(B-2) 

EAF  » 

*T 

It 

m(W+VP-UQ) 

(B-3) 

EAL  - 

PI  - R J ♦ 

X xz 

QRCI2  - Iy) 

- PQJxz 

(B-4) 

EAM  = 

Qly  + PR(IX 

- Iz)  ♦ (P2- 

R2)  J 
' xz 

(B-5) 

EAN  * 

RI  - + 

z xz 

pQdy  - ix) 

+ QRJxz 

(B-6) 

U 

■ U , u * 0 

0 

V 

" V0  * v , 

vo  - o 

w 

= W + w , 

W * U a 

O 0 0 

p 

" Po  + P ’ 

p *=  p 

0 0 

Q 

- % * q ’ 

Q0  - 0 

R 

a Ro  + r • 

Ro  - Ro 

(B-7) 

The  second-order  terms  have  been  ignored  since  the 
perturbations  are  small.  The  velocity  (U)  is  assumed  constant 
which  eliminates  the  X-force  equation  R-l.  After  using 
Taylor  series  expansions  of  the  nonlinear  terms  and  substituting 
B-7,  equations  B-2  thru  B-6  become 
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IAFy 

* m(U  B + 
o 

U r • P U a • 

0 0 0 

Ua  p) 
o o 

(B- 

■8) 

EAF, 

* miU^ai  ■ 

► P U B - U q) 

O 0 0 

(b- 

•9) 

IAL 

“ PXx 

- rJ  ♦ 

xz 

qR0n:  - y - 

q P J 
' ox: 

(B- 

10) 

5AM 

= 

y 

♦ rP  (I 
o x 

* V * P^o^x 

- I,)  ♦ 

(2PqP  - 21 

l r) J 
0 

(B- 

ID 

IAN 

■ r I , 

- pj  + 

* x z 

<ivyy  + * 

R J 

0 xz 

(B- 

■12) 

where  a*Aci',w/U  , 3=v/U  , and  W ~U  at  . 

O O 0 0 0 

Ignoring  gravity  and  combining  equations  B-8  thru  B-12 
with  their  aerodynamic  force  and  moment  equations  yields 
(Ref  5:4.113) 
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rlx  • PJx;  * ‘W1*’  * Vn'’  ’ ^blC^e 


*c„  (4-)  P * C„  (jg-)  r * C «f  . Cn  5r] 


(B- 17) 


P o 


r o 


The  stability  derivatives  and  aircraft  dynamic  terms  are  listed 


in  the  computer  program  EIGEN  in  Appendix  E. 


i 


APPENDIX  C 

Development  of  Longitudinal  and  Lateral  FCS 


Differential  equations  were  formed  from  each  of  the  transfer 
functions  in  the  longitudinal  and  lateral  FCS.  The  physical 
variables  XI  thru  X5,  Z1  thru  Z4,  Sf,  Sr,  and  6s  are  shown  in 
figures  2.1  and  2.2. 

Longitudinal  FCS  for  a < 20.4  Degrees 


1)  X, 


10(57.3' 


yields  the  differential  equation 


- 573ot-10X, 


!)  X. 


yields 


/ 57 . 3 S\ 


3)  X, 


yields 


1 5Az 


X.  - 57. 3q  - -X, 

L *- 

15U0(a-q) 


:3  ‘ TT71  a ‘ 


15U  a 


- 1 5X. 


vields 


(C-l) 


(C-2) 


(C-3) 


(C-4) 


(C- S') 


(C-6) 


4)  X4  " [0.161(X1  + 0.7(F3)X2)  + (0.2X2-0.SXS)]  (C-7) 
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-.161(3)X1-3[.161(F3)(.7)+.2]X2+3(.5)X3+X4=.161(15)X1 

♦15[.161(F3) (.7)+.2]X2-15(.S)X3-lSX4  (C-8) 

5)  X5  = F4  X1  * 2 . 5 (F3)  X4  + | (2.S)(F3)X4  (C-9) 

yields 

-F4XJ-2. S(F3)X4  ♦ X5  - 12.5(F3)X4 


6)  5 


yields 


(iOKSX  x 
l s+7TT  / A 


(C-10) 

(C-ll) 


20(KS)Xs-206s 


Longitudinal  FCS  for  a > 20.4  Degrees 


(C- 12) 


The  differential  equations  for  X X,,  X,,  and  fi  are 
the  same. 


7)  X, 


yields 


(tJtI)  (*2X2  ' ‘5V 


(C- 131 


-.6X2+1.5X3+X4  = 3X2-7.5X3-15X4 
3)  X5  - F4XX  + 2.5(F3)[X4+.SX1+.5(.7)(F3)X,] 


(C- 14 ) 


+ 2-5(F3) (|) [X4+.5X1+.5) (.7) (F3)X,] 


(C-15) 


yields 


t-2.5(.5)(F3)-F4]Xr2.5(.5)(.7)(F3):X2-2.5(F31X4+X5  - 
2 . 5 (F3) (S)(.S)X1+2.S(.S) (.7) (5) F3) 2X,+2 . 5 (F3) (S)X. 

«-  4 

(C- 16) 
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Lateral  FCS 


9)  ft  - 57.3(.12)(|«F)  p 


yields 

6f  - 57.3(2.4) KFp  -206f 

10)  zi  = = (rrT(^^Tr)T)[poa+0topl 

yields 

Z = 10  , 0 

1 5’.3°op  + TT77yP0a-l°Z 


1 


U)  Z2  ■ (kr1)  CS7.3r-Z11 


yields 


-1.5(57. 3)  r+1.5iL+Z,= 
1 *- 


121  z3  -(snr  )2; 


( ' 


yields 


-3Z2  + J3  = 15  Z2  - 15Z3 

13)  Z4  = F8[Z3  + GAIN1  Ay  + . 1 2 (57 . 3)  . 0375 


F8 (GAIN1 ) (.6)11  (tf+r) 

F8Z  , + 5 ♦ 

J 32.2 


— [ . 0375- . 65  (F7)  ] (Ppa*a0p) 


(C- 17) 

(C-1S) 

C C - 19) 

(C- 20) 

( C - 2 1 ) 

( C - 22) 

(C- 23) 

(C- 24) 

. 6 5 (F'7 1 1 pit 

(C- 25) 
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yields 


-10(F8) (GATN1) (.6) (Uo)  P - F8 (GAIN1) ( . 6) (UQ) r -F8^3 
32.2  32.2 


+ Z4  = 1.2(57.3) [.0375-.65(F7)]Po«  + 1.2(57. 3)  [.0375 


+10 (F8) (GAIN1) (.6) (U  ) 

3272- r + 10™Z3  - 10Z4 


(C- 


14)  «r 


/ 2 0 K R \ 7 
V S+TO  / L \ 


(C- 


yieltls 


6r  = 20(KR)Z4  - 20Sr 


(C- 


-.65(F7)]aop 

26) 

27) 

2 S) 
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LRTERIfiL  NODE. 
DUTCH  ROOL  ♦ 
ROLLING  nOOE  ROOTS: 


APPENDIX  D 


CLOSED-LOOP  RNO 
OPEN-LOO?  ROOTS, 
RT  RLPHfl=  2S  DEGS. 


I RERL-RXIS ) 


-0,28  -0.24  -0.20  -0.16 


-0. 12 


LONGITUDINAL  flOOE. 
8H0RT  PERIOD  ♦ 
LONG.  FCS  ROOTS: 


IRERL-RXIS) 


-1.60  -1.20  -0.60  -£.40  0.00 


0 £ 


OPEN  LOOP  ROOTS.  a - CL08EO  LOOP  ROOTS. 


Fig.  D.l  Closed  Loop  and  Open  Loop  Roots  for  Aircraft 


APPENDIX  E 


El 

- YFROOT 

Solves  roots  for  equation  2.8. 

E2 

- EIGEN 

Solves  eigenvalues  of  full  linearized  model. 

E3 

- RSVP 

Plotting  program. 

E4 

- NONLIN 

Solves  equation  3.1. 

1 


o *>  on  o nno  o o 


El 


P-’CGf  A * YFR03T  (Ivp'JT  = /8Q  , OUTPUT,!  a p £ 5 sINPJD 
ClHEf.SION  M5)  ,7U) 

PEAL  C,  1,  C, 2, 91,  P2, 37, 01,02,01,02,03 
COMPLEX  7 


*■  f<«m« 


v » * * *. 

' , * + a - » , « « » . « * » » . X ; » 4 » « « «.  , 


C TNPt'TS  TO  PROGRAM* 


U 

IX,IY,IZZ 
JX  7 

D 

P'lAX 

S 

or 
CE 
P nw 
-AS 


veld;! ty 

Mile  MTS  OF  I KEPT  I A 
prodj;t  of  inertia 
poll  = AT=- 

WAYI-J-  POLL  RATE 
SJHr40E  ap^fl 
WING  5 PA M 

chop: 

A IF  PnETSU^E 
-AST 


rr/sFr 

S.  JG-(F|--*2) 

T.  JG-(  - r -?) 

pa?/sf: 
ra:/sf: 
cr»  *2 
rr 

FT 

3_  J GS/ ( rT  * * 3) 
S.  JGS 


r STABILITY  0 FR IVA  T I V r 3 ? 


L. 

CNK 

- /RAO 

r 

CM3 

- /5-QEG 

c 

CYO 

- /RAO 

c 

C-A 

- /DEG 

r 

CMAO 

- /DEG 

» f.  * * • 4 .•  4449444^)  f l « * 4 I f * l > » « . *••  » \ 4 a.  at 

AA|«44«4^*tf444f4  44  44f4i  4 ► 4 «4  • JB  A * f ■•.  4 * 4 A f 4 4 • V 


K’0EG=R 
0=1  . r.. 

PMAX*  1 

•Js*  25.::77=38613.t  !''  = !’  S D C G.  j 7Y=>  .20  3 . : JY7  = 2-  0 . 
t L FA  = c . 3 s=  2°  3 . EO  r=  . S . r E r 1 3 . 9-  S7’  2 >•  = . G 0 J J 7 T ; =C  = -«  9 . HG 
MA  F=  C ?4  . •• 

CN~  = - . 7 * C*19  = . U C > " : *0  =-u 

C-tnr  c. *CMA=.: 33 » 

(-"ArC-A^iZ  .3 
T^a-r  *p*  -7  .3/5 

10  rf*A0.s,L  • ROW*  (!J ■•*?) 

Is  Gl=(TZ7«-i'Y-iy)  'O/(5,0,’r*«  of) 
r-2=-FE-C‘JR*P/(  u*2> 

°l  = I7?/(S-»n0AP'4r) 
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9 


n?  = -r  F-C'JR  'HJ'  ?) 

°3  = CtP“(iY-lX)  ' (P‘  • >)  / (S *OPAR*P- » 
Cl=(IY-IX«'I??)A-3/(S*r):5/,f*CE) 

02=-rp*0'io*o/(  u*  ? > 

C1  = IY/  (S*  0019*00 

C?  = -fE*  (CMO+CiAO) / ( J* ?) 

P3=  I1X-I”) * (»**  ?)  f(  <;*  pgA;?*CE)  -C<n 

A ( 1 ) =1.0 

A ( 2 ) = ( nl * 02+  0 1 ) / Ml  01) 

M3)  = (01*  Gl*91’  03*P‘»*  "‘’♦B3*  01)  / * M * 01 ) 

A (ii)  = (ClKG2*C‘**"»lt3V0?*B3*Di?>  /Hi  *01  ) 

A(F)  =<C2*0?+  33-031  / ('’I'll) 

CALL  7 PQL?  (A  » N nE 5.  7 , IF") 

IF(P.PO.O)  GO  TO  ?C 
c?INV*t"  ” 

PRINT* i ” " 

P?1N74‘,”  ** 

P 9 1 N ^ * t"  '* 

C?INT*,"  “ 

GO  TO  30 
20  no  in'*  no 

30  PRINT- ."THE  ROOT ? ARE  rOk  p=m,P,"  RAO/SEO 

°kINT*,»  •• 

PHN7-V  " 
no  i*  9 1 = 1 ♦ N0E3 
PRINT-*  ♦ 7 ( I ) 

*»0  CONTINUE 
P=F«-.002 

IF ( ° • L E • °M  Ax ) GO  TO  1 r 
1Q0  FOf’UAT  ( 1H1 ) 

FND 


t , 
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, .i  . . Mil  MMn 


C C‘  t C ('  C C l C f C C C'  c 


E2 


eROGf  AS  El *; r M : M r J r = / * C , OUT  D!J7 1 ’ s 3 1 3 s I NF  J n 
n*E  r sio'j  * < i7)  >w(i")lTfi7,i‘'»t-K(-a)*3(i7,m 

PIMEr  sign  a i hV( ir , iT)  , MKAr.rft  (i-o^j , S(  1 7,1 
DIMENSION  31  (10?)  , ?.l  ( 1 *’?) 

COMPLEX  T,N»?U 

peal  f3,F4,P',ei»e?,e?,e<« 

PEAL  Kf,^,<S,r5,t&,?7 


-»  » i‘  * M » * » • »*«  » * : *»  »«  .***«>* 

tti  **»>♦»*  ,j,  t UM  »*•*•*.•  < I f I 


r this  fpuC-ram  if  srr  up  to  do  a lin-ar  analysis,  in  rosy 
C 'XIC.  or  A P A 7 T I C (JL  * R AILORAFT  WITu  A lONGJTUOTNAL  ?Nn 
C '.ATEKIAL  FC.f.  IT  C?N  SOLVE  FCR  FOTH  -ISENV5LUcS  SNj  EIGENVECTORS  • 
p ™SL  SUPPCUT I N E5  ARr  USED  FOR  THE  ANSlYSIS. 


C "'p  THE  1 7 STATES,  THE  CIPST  V ARE  AIRCRAFT  STATES? 


r 


nEG  A 

- 3I0ESLI3  ANGLF 

- (RADIANS) 

ALPHA 

- ANGL  F OP  A IT  AC' 

- (RADIANS) 

C 

- PITTH  PAT p 

- (RAOl^NS/SEC.) 

P 

- ROLL  R AT r 

- (RADIAN5/5EC .) 

c 

- YAH  SATE 

- (RADIUS'S  EC.) 

*«E  FOLL  OWING  6 STATES  ARE  FOP  THE  LONG.  FCS 


XI 

3H  Y S.  CAL 

V A ?I A 3L  E 

X2 

PHY  SI  CAL 

varta^le 

X 3 

°hy S ICAL 

VARIABLE 

XA 

p h y s TC  A L 

VARIABLE 

XE 

fhvsiCAL 

VARIABLE 

0E*.T  A -S 

- GTuPlL 

;ppR  AiGLE  - (P 

0 

” H E L **ST  o 

5 7 A T ES  A 

r E FOR  THE  LAT. 

r 

DEL  T A -F 

- F L fioER0N  ANGIE  - (PAD 

e * 

71 

PHYSICAL 

VARIABLE 

r 

7? 

°HY  $ r A L 

7 A ?:  A °LE 

•?  7 . 

PHYSICAL 

7 A R I A n L E 

+% 

7,» 

PHYSICAL 

V A R I A nL  E 

r 

DPLTi-R 

- RUDDER 

ANGLE  - ( RAO 

«»• 

^14444-W'  ft  f . * * * : * + «*»+*•+* 

L66 

...  i— qTcgrojr— T— «t_ 


c 


C "’JPtjTi  •'0  P-.;OOA  AM! 


r 

U 

- 

vplo:i ty 

- 

r V / S EC 

C 

alt 

- 

alti r jor 

- 

FT 

r 

IX, IV, 127 

- 

MOSERS  OF  INERTIA 

- 

S- JG-(CT* “ 2) 

r. 

JX7 

- 

pROO'ICT  OF  INERTIA 

- 

S _ U G - ( r T ' *2) 

r 

ALFA 

- 

AN'GL " OP  ATTACK 

- 

O'.CRFrs 

C 

P 

- 

ROLL  R AT  r 

- 

RAD/SFC 

c 

PMAX 

- 

waXT*^JN  ROLL  RATE 

- 

RAO/SEC 

r* 

P 

- 

YAW  -ATE 

- 

RAO/SEC 

c 

S 

- 

SURF ! 0£  APPft 

- 

FT*  *2 

c 

CE 

- 

CMOR  0 

- 

"T 

c 

nr 

- 

UNO  SDAN 

- 

- r 

POP 

- 

A I P PRESSURE 

- 

S.UGS/ (=■!«»  3) 

/N 

l , 

PC 

- 

STATIC  3RF^SURF 

- 

uBS/  (FT'  > 2) 

r 

HAS 

- 

MASS 

- 

S- JGS 

r. 

~TABItn  Y 

DERIVATIVES  ! 

C 

CNC 

- /RAD 

CYR 

- 

/RAQ 

r. 

CNF 

- /5‘D PG 

CYB 

- 

/DEG 

C 

CNP 

- /RAO 

CYP 

- 

/RAD 

r 

CNDF 

- /OEG 

C Y OF 

- 

/ 0 F G 

r 

CNCR 

- /DEO 

C Y DR 

- 

/ 0 E G 

r 

CLP 

- /RAD 

C7  A 

- 

/RA  D 

r 

OLE 

- / 5 1 D EG 

C7AD 

- 

/RAD 

r* 

CLf 

- /rad 

C20 

- 

/RAD 

r> 

Cl  CF 

- /DEG 

C7DS 

- 

/RA  D 

c 

CLDC 

- /DEG 

CMA 

- 

/ 0 FG 

c 

CMC 

- /RAD 

C^AD 

- 

/DEG 

CMCS 

- /DEC 

r 

P7,F/  , 

F7, Fl  - -CS 

nAr  AMET  EES , 

s 

0.7ED  r OR  7 N PROGRAM. 

C *F,Kf  tK'S, 

c GA1N1,GaTS?  - Dll  Ml  Y CARAME7C RD  JSED  TC  A V A L Y S E CERTAIN 

C ““"I A C < LOOPS,  J“’I\LLY  SET  TO  1. 

J l V ■ » » » .*  •.»  r «.  .«  < * * » > - > '?»*  > • + *-  4 * » ••  * • 4 '*.»  4 • , i 

« •"«  « *>•#.>  K t * * « a +-•*>*.  I - R *.  * ».  «M*4,  ».  ■»,*•»♦*»  I**'-  ~ 4 ■ ; •*-»"■  * *•  ' • 

C The  program  runs  for  the  last  value  of  alfa  read.  Interchange 

C cards  to  run  for  different  values  of  alfa. 

ALrA ~ 0 • 

A L f A = 10. 

row=  . C 023-3  I PC  = 211  =,*  ALT  = 0 . 

. SI  ?7  = SrfcOO. : IY  = l"C00.JIX  = -.'>r  J.  i JX7  = 2«-  G.  S&L  T = 3M.0D. 

ALFA=?C.  -S=  ?57.>:E  =10  . ^ 4-Br  = 20 . • RD*l  = . 0(10  7 V 5 = C = :+p?  ,qr. 

C-PAr?*.  . 
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VAS  = i 3*-. 

ALFH* 

1 Ic ( ALC6  .ST . 0 ) GT  T T 7 
IF (ALFA .£">.  1C)  ZD  TT  3 
IF(ALFA.ET.O'i)  ZD  TD  *- 

CNR=-.‘,0VK  V*=.C  VO- r>  = 2y*C'1*  0=.  29>2 

CMA  = . C Cf  i CL  n F=  , I 03  3?.'C "0?=-  . 0 C 0 1 f OL  TR  = . 0 j 0 : 2 
CNOR  = -.OiU  j?  :TSD:.0'’48<*CLP=-. 13?^  00: 

CLP=  .271.-*  ?cv  P=-.  r l F ♦ ?cvne  = . 00001*;rT\  = .C030S 
r7A  = -3.^ifT7/;-'  = -l.  5177-C?0=-3.2^3:  :TDS=-.  S767 
ni0S=“.0G75?CYo=0«  STYP^o, 

r-o  t r 5 

’ CNi\=-  • if  3; : n = .oir*  rT^n--^ i c^ia'y = c . 

F^As.GPCjSCL  C”=,  r2lT-j?rNDF=.  OOO?*OlT9=.C3057 
CMC.R=  - • PC  1 ■*  i r;^=C  . !3LP  = -.2:SCL',=-.  00  4b 
CLR=.  128iCvP  =-.0  23^.3':  Ynpr-.  c 002-  : :YT?  = .0  0 334 
C7A  = -3.9;C’AT  - -i  ,(?o  VC7'->  = -2.0G7':733  = -.?  73 
C*iTS  = -.  OiGOY^Q. ;OY°=P. 

cn  tc  s 

7 rMks-.KOUiCf  3=.ri3>TM,'=-c.7C0F5:-1.\3=.lC0? 
ci  a=  . fj  or.7  ? : l:-c= . : o 1 1 scnd^t  . oo  c * > i z _ tr=  . c o c b-i 
CM9Ps-.0Gi’r*3N»s.  au  sn.°=-.2se*  * 3=  - .c  077 

CLR=.  227.  3YR=-.0  2?33CYnF=-,oorir3Yrl^s  . Of  7?^ 
r7A=-3.C"o95C73T  = -l.?F£bfC7T  = -2.;tt3i C7P5  = - .bV 
Ci3S= -.31*0  >&=0. ■ 7yo=C. 

r,o  tc  b 

n NS*  - . F»\d  * * c kir*= . 0 0 3 ? * co  = . 2 1 «,  9 ■ r <i  n = - . 21 4 3 

CMA  = -.  0 03  5?L3C=.  JO  0 ♦ 2?rN0Fr-  .C0"l+i3L  HR* . 0 J P-.  9 
CMCP  = -.0C136'CN=  = . n37:'CLP=-.l?72r:.c'  = -.C0?5 
CLF=.?3l7  3CY  ==-.  r 1 * l .*  C v PF=- . OOJ?S:TYnR  = .j0?d.* 
C7A  = -7.O72:C7AT  = -?.l'>g-;C7''=-H.',:33J3  705  =-.*>73 
CMCS--  .OC  7 i£CYP=; . ?TYP=0. 

' ITGT  = 0 

row  COfiV-Rr  ST  A 3 • DERI'/.  TO  TESTPcf)  JV1XT>. 

C7r S=C?OS/  -.7 
C-s)f»=f  K.‘",  7.  VS 
CL°=f  l?-, 7.  7/r> 

CYr,=  cyn»;  7.3 
C*ia  = f . 7 

iJOa=C.  P'JP  -IS-^7AL('r5  CMLY. 

IJOq*i.  pop  £T  "*E  M')A  L'lSS  A»jn  FI  * r c Tr  R3  . 

I J0«=1 
? jnn=o 
c=u 

«-‘*AY=12. 
p=c 
: 1 = 1 


1 


*=!7 
1=  17 
17 
N=  1 7 
TA  = 1* 

13  = 17 
IC  = 1~ 

17  = 17 
KF  = -1 . 

KR=1  . 

*S  = 1 . 

GAIN 1=1, 

GAIN  7= 1 . 

ALFAr  At  F\/}7 
GlAR=.'  * -ON-  ( l • * - ?) 

F ATI  T =n$UR/ 

CAtCl'l  ATjONS  F07  FCS  P A RAMFT  EkS  • 

Fl=(-.6V  (73AR-1LO  + .3C. 

r?=(  ?2**<  n$4*v-'  o :>  /77  0 J+.22A)*.  1 

F3=  .*  ♦ ( RA  T I 0 -7 . 3 <» ) / 1 . 1 ° ? 

FH=t-l.*?M:?ATI0-.57)/i.?t,*l,F)-l.n 
E<5=1  . Ct  1«-  MATI0-.1P*? 

F6=- ,f  F7*>»-  <A  1T0*?.127fc 
IF ( 0 F A R » L T • 1 D *3  = 1,  C 
IF  (PfAP.GE.  1 r0>  r7  = -71 
IF  (O^AA  ,Gf.  C-0  !l)  r-*="2 
IF(QrAP.3F..3:o3)  f3=,176 
1F(p7TIG.LT  •.rj)  Ci~.T 
JF  <Fp TIO.GF.  .7  ?)  -A:-4 
IF(Rf  T lO.'r: . 1.7  31  F*  = -1.0 
IF(PAT  IO.LT  . u.IA;  ) <-7=0. 

IF  (RATIO. GI.  P.  H.  ) r,  = c(> 

I F ( F.  f 7 1 0 » G ” • 1 . 1 7 ; ) "7=1, 

Tr  (FAT  xO.GF.  1 . □ t ) -7  = r6 

Ic  (FATIO.r,-:.  T.  ?■»)  r'  = t>. 

:c  (F  OTIO.LT  .2.  CU ) F-,  = .“ 

]F  (pf  Tio.i*  .?  . 3V>  F^  = r-» 

Ir (RrTTO.GF. *.23)  nn,0 


5F1  t*P  STWr  ^4fvIV  f A)  *XOOT=(P)  • X • 

A*F  tJL  V"  r r < El  C,rN  VALUES  PF  (A-TMFP.SO  »3 

♦ 

♦ 

0-»  #••••»>  * * * • * 
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n0  1‘  1*1,  *1 
on  i»  j=i*n 
a (i , J) =o 
° < r , j)  =o 
1*  CONTINUE 

* (1,  1)=MAS'  II 

A(2,<)=*A  S'  'J-0’>4r*S*C’!\0 

t ( ;,  ?)  = -rruR*$*  ( ;■•*?)*  cmao/  (?»'n 

A f 5,  7) *IV 
A (*  , * ) - IT 

a <-  ,*  )=-jx7 
A (C,t.  ) s-JX7 
A (‘j  ) = I77 

no  2'. 

A(i*I)=l 
?’  CONTINUE 

A ( 7 , 3 ) = -•>  ? . 3 

A ({?,  /3Z.? 

A (?,» ) *-.*i  3 

A(9,?)  = (.tol  - .7*  r 3 t . ? ) 1 ( - 3) 

A(<r,f  ) *1.3 
A l 1 C » 6 ) =-c-» 

A.  (10,«>)=-2.f  «F.T 
AC  J 

a(1h,i7:=i,e 
Afl»>,10=-3 

e ( it, i)*-n»PH i Ni*  .r  ♦u/.^z.? 
A(16,5)*-c9*3AINl*,4*U/32,2 
A ( 16  , ID  *-rt 
P(l,  1 » - 0‘3  A-R  * S-CY3 
« { 1 , ?)  =MAS*  U*  P 

0(1  ,M  =N*.  3*U4  ALe6  OHIV  S*OF*CYP'  ( ?*  J ) 
P<1,‘  )=-*n  ;*lH-OtU*'V  oc*CYr/ <?*U> 
o(l,l?)*0aAO-S-*3Y0P 
P(  1,  U ) s^K'  S*:vT--, 

r rr, 

0(2,  ?)  *Or'l 

r,(I,3)  = l*a^i'  + 0"»:?'f,‘Zr  C ?0/  ( 2*11) 

O < r,  1 1 ) *nyU  *9-  0 DS 
r(  j,  ?)*-OrA7*5-or’  Z"* 

° (3,  ’>  =01 A <'  S»  C»  • »>  * CM0/(  ?*U) 

*•?,.  )S-C«  ( IX-I7,»J  -7*  jy7,D 

">(3,'  )=,'”■  <- jy^-a-  ( I x - 1 ** 7 ) 

r(3fll)sO”Ar»S'*“r,C1T^ 

0 (‘  , 1 ) =0'D  * lc  I L T 

P(.  , U.-P  JX7  -0  (’  77-rv> 

n(,  f , n“*'?)*"lP/(2*U) 

~M  ) =1*11  ;♦  .;•*  r*c«  • ■»)  ♦ ^LP/(2*’J) 

° U , 1 2)  *00AC  -S  C-l* 

«(•.,  lt)=n-nr‘  <>■)«■•  01. op 
n (:  , t ) spo.lv  ► s*  *MT 


-0 


?<!,,  ?>=-P*C  >-i  *>  - R ' JV7 
P(F,<  )=QPA\»S*  (IF*  *2)*''NP/(2*U) 
r<£,!  )=Or«A^*SJ‘(1F»*2)  *rNR/<2*U) 
p ( F , 1 ? ' =fi P A K • S’  3 - » C 'IT P 
P(i  , 16  ) = C3^  R * 3*  Pr*CM3P 
9(1  , 2 ) = 10*  57  .’ 
n<fi,f  >=-10 
PC  ) = -l 
B(8, 3>=-'J*lF/72.2 
,f  >=-13 
«<9,C  >= 2.415 

P(9,  : ) = (.  lb  l-.r*r7*.  2)  15 
9<9,f )=-7.0 


ft  ( <3  , « 

8(1G  * 

R ( 1 1 , 

R(ll, 

«(12, 

9(12, 

9(13, 

9(13, 

9( 13, 13) =-l 0 
9 ( 14 , 1 4) =-l 
R ( 15  » 1 '•)  =13 
9 (1*3,  IF)  =-15 

°( 16, 2) =1.2*5 7.3 1 ( .r»7-.-.6,  * R7  > * 3 
n ( It  ,4)  = 1. 2**57  .3*  ( .T*7'  -.b5*F7)r  l.rA 
»(lo,F)  = 10‘F  ^OAm^.ft-  U/32.2 
r ( 16 , If >=13* F8 
p ( 16 , 1*. ) =-lC 
9(17,10  = 20* 

R ( 1 7 , 1 7 ) =-  2 C 

rFCALPA.Lt.2C.**)  30  TO  28 
A (6,M  = C. 
t (9,7  > =-.3 

r ( to ,l>=-2. : ■ c7-- •-*■ 
ft  (lu,7>*-2.;.  • . i*  .7*  ( r 3“  * 2) 

p (<r,i  » =0. 

9(6,7>=3. 

9(l0,f  )=2.  5-5*  .-••  7 
« ( 1 C )=2.b*  .•-•./•  5 » ( ct*.  ?) 

call  i:fme  ( : , s , i a fAiNvfior>T  » wo  l er) 
oo  3 c :=i,  ( 

"'0  ?(  J=1,M 
ft  (I,  j)=Aifj>;  ( i,jj 
reM'NUF 

TALL  V^ULP-  (A,P,.,'*,'JfTA,I'>,C,l'“,rFR) 

?0  37  1=1 , \ 
p0  37  J=1 , M 
A 'i, J)*C(I, J> 

CONTINUE 


) = - 15 

9 ) = 12 . 5 * f 3 
1 G > = 2 J * <S 
11)  =-20 
O =t  7. 3'  2 .* 

1?)  =-2  0 

2) =,7^*3 

4 > = 10'  5*  .7*ALFA 


7 7 7 ' . 6 . * E7  ) * 3 

7 77 ' - . br  * F7 ) r A _ r A 
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CALL  F1G*F(  a ,S,n  ,1  )0P,W,7,I7,Vi*,IER) 

IF<P.EC.J>  rO  TO  T5 
PRINT*,"  " 

PP 1 NT  * , " " 

PRINT*  ** 

PRINT*,"  " 

GO  TC  i-5 

FRINT  100 

PRINT*  , »••»***•#»•**•*♦*»♦♦*  *lOOr.  • +•*»-»  ♦ 

PRINT- ,"  " 

PRINT-,"  " 

PRINT- ,"  " 

PRINT* , "FOR  VELOri  TV,  N»",'J,"  Pf./SFC." 

PRINT*,"  " 

PRIfJ*  * ,”AT  £ N AL  TI  TJO  F OP  ".ALT,1*  =7." 

PR  I N T * ,"  " 

PR  I N” * , "AND  AN  Or  ATTACK  0-  " , A L Pu  A , DEIRrEE" 

PRINT  * ,"  " 

PRINT* " 

PRIM’" , "PEE ORACK  3 ARAMCTE?> « ” 

PRINT*1  ,"  " 

rG ] NT  * , " F ?=",  = ? 

PRINT*,"  FA=",PV 
PRINT*',"  Po=" » F5 
PR  IN” * ,"  " 

PRINT*  •• 

PRINT* ,"  " 

PRINT* , "THE  EIGENVALUES  ARE*" 

PRINT* ,"  » 

PRINT*  , "FOR  <OLl  3ATt,  Ps",P,"  VO'SEP" 

PRINT  * ,"ANO  FOR,  TAN  RATE,  RAO/PEO" 

PRINT*,"  " 

PRI N’ * ,"  " 

00  p>-  1 = 1, N 

PRINT*  ,W(I) 
rCNTIHME 

IF<I JOP.PO.P)  GO  TO  ? C 

PRINT* , " » 

PR IN’* ,"  " 

PRINT*,"  " 

PRINT*  , " f H E 0 0 ” 5 P T N D T N G El  GEN'* : I r D R R £R't" 

PRINT*  ,"  " 

PRINT* ,"  " 

JJ  = 1 

JJJ= JJ*3 

oo  hr  1=1,  n 

WRIT  c ?0C  , ( 7 (I,  J*  , J=  J J , J JJ) 

CO*  T T SUE 
PRINT*  ,"  " 

PRIN’* ,"  " 

J J=  J JJ*1 


■ - mm 
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ip < jj.lt .m>  go  r:  it 

CO  6i  1=1, M 
WRITF  ?2D,?(I,N) 

: CONTINUE 

*S=1-KS 
K*=1-KV 
Ke  = 1 ♦ KF 

TF(KP.cO.O)  CO  TO  o 
GO  TC  302 
F = P».2 

PsF^TANtttLPfl  ) 

ic (o  .LE.PHAX  ) G3  TO  : 

1"P  FOPNAT  ( 1H 1 ) 

?nr  FOPvrT  (IX,  . (1  M(,  't  ?.  =,  , ?H,  ,Pl2,r  , 1H)  , 5X)  ) 
?-"*  FOPMfT(lX,lH  l,Gl?.6*?H»  ,£12.o,H)) 

30.'  STOP 
FNO 
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E3 


PPOFF AM  (I\,3M'P,3'ITP'JT,TAPE71 

PlPEfT.ION  A(b2),?(£'»),P(^?),0(5‘,),'(P  2),r(r-2),5(=2),H(52) 
PIwEf-SIC*'  PI  (5  ?),  °2(r?)  , P3(b2),P  + (5?)  , Rif  i2>  , R2(  C 2)  , fr  3<  j2) 
DIMENSION  A A (104)  » H 1 ( i (?*♦ ) »PC  ( 1 O'  ) , 01  ( 1 0*  ) 

DIMENSION  tEUQM  , ecl  10<.)  ,r.G  (10-  ),HH(1CU 
DATA  N,XLEN, YLEN/; D,?.r,2.E/ 

PATA  V,V1,'/2,VV7.,1.4,t.L,4.8/ 

PATA  AL1,AL  2, Al 3,414/0 . ,10., 20.,  ?5./ 

PATA  A/50'1.  ■!/ 

ALl=/Ll/t7. 3 
AL2=AL2/^’.3 
A13=AL3/F7. ? 

AL^sAL-/1;/.  7 
R ( 1 ) — C • 
r (i) =2. 

E(l) =2.4 

r,  ( 1 ) r 2 . 2 

C1  (1)  = 0. 

=2(1) =0. 

°3 ( 1 ) = 0 • 

P4  ( 1 ) = 0 . 
no  1C  1=1,7 
B(IM)=9(I>  ♦ V/"Q 
C(IM)=C(I>  4V1/5C 
rm=7.*c(i>/i.'f-m 
E ( 1+ 1 ) =E( I ) +V2 /5T 

TF  (E  (I  ) ,LE.«  .<»)  F<  n =3J  (E(I>  - 2.n  /2 
IF(Ef  I)  ,C,T.L  .4)  P(I)=E(I)-1.4 

ZF(E  ri).GT.(  .4>  C(D=F. 

R(I  + 1)  =G(I)  + V3/5" 

IF  (G  (I)  .LT  . 7.4  ) H(I)  s6(I ) / l.  2-2.  2/1.  2 
IP  (G  fl ) .GE.  3.4)  J(I)si. 

Pl(lM)sPKI)  +7./?n 
F1(I>  = P1(I) * TAN ( AL  1 ) 
r2  ( I ♦ 1 ) ="J2  ( I )♦-  ./il 
F2(I)=rC(I>  *TU(fi?) 

D3(i+i)==3(i)  f/./n 
P3(!'  = ?7(I)'  TAN  ( f L T) 

(I*l)=°-  ( I)  «•/./;  A 
F4(I) = °4  ( I ) * T AN( A.  4) 

1*  COMTiMJE 
no  2r  1=1, N 
AA (1 ) = A (T ) 

AA  (I  +ED  -31  (I  ) 
a?(I)=°(I) 
nB  (I  *‘<i?)=Rl  ( I) 

CC  II  ) - C (I ) 
rc  ( I *C C) =°2  (I ) 
ni(i)=nc) 

PP  <I*t  C)=f??(  I) 

F E ( I ) = E ( I ) 


j 
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FS(I*5C)  = °T  m 
CF ( I ) =F ( I ) 

CF ! I *5  0)  = R3  ! I ) 

0G(X)*G(1) 

GG(I*SO>*P*»  (I) 

HH(I)=H(I) 

HH ( I *5  0) = R4 (I) 

' 'n  CONTINUE 

CALL  PLOT (?.  ,-2. , -T> 

CALL  PLOT  < T.  ,1.5, -3) 

CALL  PLOT  ( 3 . ,6.5,-.?) 

CALL  PLOT ( b . 25,0. ,-?) 

CALL  PLOT!G  -“») 

CALL  PLOT(-e. 

CALL  SYM?OL(  2,'l,5.,rc,.03,?lHa''3T  Vi  Vr_3:iTY  AXIS, 0., 2D 
CALL  SYMBOL !2 .09,5 ,?r , .36, 22H-  LIVES  OP  INCTA?IL1TY,0.,2?) 

CALL  S > V P OL  (2.??l,e.12‘ ,.05,30H-  LIMES  FOR  VELOCITY  AXIS  ROLL, C., 3 
♦ P) 

CALL  S v Ma OL  ( l.??r,3,t¥*»5,  .0  3 33 , 1 1,0.  , -1) 

CALL  SYMOOL  ( 2.  D1  , e.16  *> , .06  33, 11, 0 . , -1) 

CALL  SYVqx  ( 2.37r  , * . C , . 0 3 , 1 9 H A0  A r A V', L E Cr  ATTACK, C. ,19) 

CALL  PLOT  ( 1 . CP’*,  ► . 219, -3) 

CALL  PLOT (tpf 3., -2) 

PALL  PLOT  ( - - . 1 9C  , -3 ) 

CALL  PLOT ( .K  ,0., -2) 

CALL  PLOT  ( - 2 .2  *-.0  , -3 • lh-  , -3) 

CALL  PLOT(.«?r., .o,-?) 

CALL  SCALE!  EE  ,XtFV, 100,1) 

CALL  SCAlc(cP,YLEM, 100,1) 

E (N+  1 ) =EE  <1  0 1) 

E ( N ♦ 2 ) =EF  ! 1G  2) 

P3!N+1)=EE! 1 01) 

P3 (M *?) =EE( 1 02) 

F <NM)=F-  (10  1) 

F ( N * r ) = FC  < 1 0 2) 

<?-»<NM)=Ff(1  Cl) 

P3  !M  *2)  =rc  ( 10  ?) 

CALL  axis  (i) 3 .,  1 OH*  (OAO/SEC)  ,-10,  X.EN,C  . , EF<  131)  , EE  (102) ) 

CALL  AXIS(0.  ,0.,10HR(OfO/SFC),ln,YL“M,9C.,FF(l'Tl)tFP(132)) 

CALL  SYMBOL  ( .1?K,  9 .5,  . , 12M20  OE3.  AOA?,0.,1?) 

TALL  L INE  (1.  F,n|,  i , 0,  0) 

CALL  LINE(®3  ,R3,M, 1,10,11) 

TALL  PLOT  (C 

CALL  SCALE! Al ,XLrM, 100, 1) 

CALL  SCALE! 9a,VL“M, 130,1) 

A !))♦! ) = AA  ( 1 01) 

A ( N ♦ ? ) = AA  ! 1C  2) 

PI  !N  * 1 ) = A A ! 1 Cl) 

PI (N ♦ ? ) * A A ( 102) 

« ( N ♦ 1 ) =RT  ! 1 0 1 ) 

P !N»D  =90  ! 1C  ?) 
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PI  (N  *1)  = n ( 1 "1> 

PI  (N  + 2)  =V5  ( 1 3?) 

CALL  AXIS  ( 0 . , 0 . , 10M->  (r?T/SrC)  ,-10tXwCN,C  .,AA(101),AA(132>) 
CALL  AXIS  (1  . ,9  .,  12  HR 'FAD /SEC)  , 10 , r L E N , OP  . , 3 1 ( 10  1 ) , 3*  ( 1 02  ) ) 
CALL  SYM^jL  ( .1  JG,  2 • c » . r * » i?H  0 TPS.  AOAi,0.,12> 

CALL  LlNEd,  E,N, 1,0,0) 

CALL  LINF  (3i,  R1 , 'J , 1,  10,11) 

CALL  FLCT(3  . ,0  ., -X) 

CALL  SCALE(  :C,XLrS, 100,1) 

CALL  S^ALFt  10,  YLFN,  100,1) 

C (N+l )=C0 (101) 

C (N*  ?) =CC ( 1 0 t ) 
r2 ( N + 1 ) =CC (1(1) 

P2 (N  *2) =CC ( 1 0?) 

P(M«-i)ann(io  i) 

D(N+?) =01 ( 102  ) 

R2  (N*l  > = 00  (101) 

R2 (N*2) =90( 1 22) 

call  AXIS (2.  ,0 10«3(ean/S^C) ,-10,X_EN,C.,QC(101) ,CC(132)) 
CALL  AXIS ( 9 ., 0., inH^(DAO/SEC) , 10,  YLEN ,HC. ,20( 101)  ,00( 132) > 
CALL  SYf'OOL  ( .1  2"  , 2 .P  , . PC  , 12H10  0 F3  • AOA»,0.,12) 

CALL  LINE (C ,J , N, f , 0, 0) 

CALL  LIKE  (32,R?,'i,  1,  19,11) 

CALL  PLOT  ( 0 -’) 

CALL  SCAlF(GG,XLEN, 130,1) 

CALL  SCALFMh,  YLE^,  123,  1) 

G(N>1) =GG( 101) 

C(N*2)=GG(1C<> 

P4  (N*l)  =C,G(  1 Cl) 

FA (N+2) = 3G( 1 02) 

F*(NM)  =HH(  10  1) 

H ( N ♦ ? ) =HH (102) 

R4(N+1)=HH( 101) 

F-  (N*2)=-H(l  02) 

CALL  AXIS  (0  . ,0  .,  10HP  (P'C/SPC)  ,-13,  Y.PNjC  .,GG(1  j1)  ,GG(1  32)  ) 
CALL  AXIS  ( 0 . ,0  .,  lOM'?  (R-'O/SFC)  , 11,  *LEN  ,l3C  . ,HH  ( 101 ) , HH  ( 132)  ) 
CALL  SYMBOL  (.  1 ?S  , 2 .•■>,.  n , l?H  2*>  0E3.  eOA*,C.,12) 

CALL  LINE (GjP,^, 1,0,0) 

CALL  LINr(PL  1,10, 11) 

CALL  PLUTE 
FNO 
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P?COf  A*  NO‘Jl  lNC\oiirtoi!TPUT,Tfl°c’rtO  JTPUT) 

C0*MCN/CALCUN/C1, C? 

DIMENSION  < (?)  te  « *») , A JTNV  <?, ?)  , M < 1 3 0 ) , Pi  < 3 0 » 1 00  > » S C*  0 , 1 0 0) 
DIMENSION  AA<  1C  3)  , =»3  ( ICO)  ,CC  <lOr>  ,T)(  IOC  ) , EE<  100) 
nI MENS  ION  FF  (130)  , GGM  CO)  ,UH  (10t)  , COUNT  (30) 

OIME*SICN  OX  (1) , XMAX ( 1 ) , XMIN ( 1) , INI T < 30) 

real  init,ri 

IMTEf  FP  Y,V  ,11,  Jt,  VI,  V5,  V3,  VA 

1 * » “ . ' • «««**  ♦ +■***..*  » 4 .1-  * * * * * * »**•••#*««  ♦ » *■♦*♦♦*  **  »«"'-».*'4'*'*»**-#*'*-»  «•  » * .»  f « » * 


0 

INPUT; 

TO  PPQGRAV 

r* 

Cl,C? 

- STARfjNr.  v£LUE  Or  CONTROL  VARIABLES. 

C 

1 1 , I ? 

- INTEGER  COUNTS  OP  INC?  EMENTG  0F  Cl  E C2. 

n 

N 

- NUMBER  °F  EQUATIONS. 

<** 

D,P 

- INCREMENTS  OF  Cl  f C?t  PFSPCCTI\/EL  Y . 

- 

F»CP 

- THE  * OP  - VALUE  TO  ATC^T  FO?  'tIFUPCATION  POINT. 

c 

INIT 

- HE  INT”IAL  G'lESP  X PC  R EACH  INCREMENT  0C  Cl. 

c 

ClMAy , 

C2MAX 

- The  MAXTM'tM  VAL!IC  OF  Cl  r C?. 

r 

PSTEf 

- STE0  L F‘‘GTH  USED  IN'  NSOIA. 

r 

M4  XF'UJ 

- MfXIH'IM  NUMaEP  Oc  EXECUTIONS. 

P-AX 

- CST IM A T c OF  r* ISTA N C r TO  C-'LU7I0N. 

ACC 

- A-C’RACY  PEQ'JIRFO  cor  SOLUTION. 

* ¥ * 

i4> 

M 4 4 4 + A # -C  Jl  # * 4 * 

PAT  A n,C?,  1 1 , J1  ,»ltTt  0,  OFLTA/-^.  , 1,1 , ?,  .r',  .2,  . 001  ' 

PA  "i  A 1N1T/1  .A,  1.  .*  , 1.  ?.  , l.R , 1.  <■  , 1.  7, 1.  ? , 1.  1 , 1 . , . 01  ,C  . / 

TATA  ClMAX»C?MAX,X/X.0»,.»?.»2./ 

DA*  A OST  E3, MJ  vc  JM  f T(AX,ACC/.001,10J,r0.,.lr-l0/ 
r A TA  FcOk/.(iE/ 
p^iN-  ic: 

FR1NT#, *»**'■  **OAT‘  FOh  SINGLE  NONLINEAR  ETJf  T TON1''*'  '*  * 

PR IN’* ” 

P’lN’’  * t **  ” 
p.rIn  • « ,*•  •* 

C-tLL  SSOla  (*<  |X,ct  A JTNVtOSTCPfOMAX»  ACC  XFJ?4, 0»  W) 

*l(Ii*Jl)*X(  1) 

tall  pknT(X,  e,-o,fi,r?,N) 

:r<ci  .'*E.o.  o go  ro  10 

r ALL  rELI^(X,Cl,3rLTA,vMAX,xy.IN,Fl,E''0o,*U*?IFr*c»Il»  JlfXA  t) 

1*  ?(I?.»J1)*C? 

r?=c?»P 

iF(cr.  "T.p’hix)  :-j  ro  *o 
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4- 


.......  „ 
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JlsJl+1 
r-o  to  • 
ci=ci+o 
02=-’. 

COUHT(Il) = Jl 

IF  (C!  .GT  .C1MAX)  r-0  TD  r5 
Jl  = l 
T1=I 1+1 
X ( 1 ) -INIT(I1  ) 

X ( 2 ) =INIT(I1) 

GO  TO  b 
• STOP 

? CALL  FLCTT(P.  l>S*AA>Tr5ff‘C»OD»EE»!rr»S3jHN,33Jf4TtT  j) 
r FORMAT  (1H1) 

END 


SUBROUTINE  CALF J^(M,X, E) 
C0MMCN/CALpUN/Ct»3? 

Cl. YE  NS  I ON  X ( 1)  ,=-(1  ) 
F(1)=02MT(  1 ) * X ( i ) +T1V  X ( i) 
F (2)  =C2+  (X  ( 2 )*v  ( 2)  +T1)  ' X ( 2) 
PETUf  N 
FNO 


FSO=FS0*F(Y)  ‘P(Y1 
1 CONTINUE 
PRINT*, “ " 

PRINT*, *'  " 

PR  I NT  4 | *'  " 

PRINT*,  ••  •• 

PRINT* , “FOR  ri=“,:i,”  1 C2=“,F? 

PRINT  ?00,  < V,X(Y>  ,F{  Y)  , Y-=1,N) 

PRINT  on, 

?1  " FOKMiT  (//^X  ,“Y”,  ! v ,“Y  ( t ) “,12X,  “r(  r>  ”/  / (I  j , ?F17.  T ) ) 

?1"  PPPM.'T  (/i-X,  “TM  E SJM  OF  SQUARES  :>",E17.f) 

RSTUr  N 
cNn 


f U*TF- CUT  IMC  P (X,  11  , ttlTA  , Y^AY,  X^IN,  PI  , r:OR,TY,niPF,?,l , .1,  YAP) 
nT“PKSlCN  Y ( ?)  ,Y“AY(t)  ,XMI‘  (l),*l(T0,10n,TX(l) 

py ( i j * c i / ?)  *'.* 

Y MftV ( 1 ♦ 

YRiN  ( 1 > =T X ( 1) -£R3R 
ya pa; p${ y ( m 

ic  (X'  a .LF.  YMN.O'  . VAR  . V»14X)  * * T J R V 

PIPF-F.l  (I,  J)  ■OKI,  J-l) 
y ( i)  =x  ( i)  «.-)tff-.  . • iTfr 

X ( 2 ) = Y ( 1 ) 

Ps -P 

PETUf  N 
PNO 


-0 


I 


S’JflRC'UT  IiJE  P LOT’  •?  1 , 5 , A A , OP  , CC  ,n 3 , EE  , Ec , r 0 , HH , 3 0 L'N T , I f J ) 

PIFEnSIOn  ^1  (’0,1QP)  ,S(30,10H)  ,iWlDC)  ,n(io?)*cc(ico 

DIMENSION  00<  130) , EE<10  0> ,FF (131) • 33  < 13C  > ,MH(  130) , COUNT ( 30) 

IN7EGE*  VI,  V2  ,73,  VR 

XLEN-2.E 

YlEN=2.S 

Vl=COUNT(l) 

V?=CCUnT (3 ) 

V3=CCUUT  ( 9) 

V4*CCUNT(13> 

PO  3 0 0 J=l,  VI 
SA(J)=F1(1,  J) 

PO ( J) =S ( 1 , J ) 

CONTINUE 
r 0 3 ? 0 J=  1,  V 2 
rc ( J ) =R 1 (5 , J ) 
co  ( j)  =s  (> , j > 

CONTINUE 
PO  320  J=1,V? 

EE  (J)  =c.l  (9,  J) 

FF< J) =S(9, J) 

CONTINUE 
PO  330  J=  1,  V** 

GG  ( J ) =R  1 < 1 3 , J) 

MH( J) =S( 1 3, J) 

CONTINUE 

CALL  eL0T(2.,-?.,-3) 

CALL  PLOT ( 0 • , 2 -3  ) 

CALL  DLOT<- .5,-.?; ,-^) 

CALL  PLO” ( 0 • , 6 »7 : , 

CALL  PLOT  (o  . 21*  ,0.,-2> 

CALL  PLPT  { 3 ,,-J. 

CALL  PLOT(-fr.??,C.  ,-?) 

CALL  FLOK.  3 ~,.*0,-7) 

CALL  SYMBOL  ( 2 • ■"  • **  1 ■* , • OS  , 16PVI D T = X-*  3*  ~V  X«-Cr,  3 . , 11) 

CALL  SYMBOL  (2 . 12?,  . 3?  , 19H-  BIFURCATION  NT  , 0 . , 1 3) 

PALL  Sv-->Ql  (1.  i ? * F,’.1  'o,.Cd,11.0.,-J  ) 

CALL  SY’^OL  ( .J  *2'  t 2.y  , 1uy,  3 . , l> 

CALL  SV;'jOl  , .lsb2,lHX,  3.,  1) 

PALL  SYMBOL  ( 3.  3‘> 1TC2,  1HX,  e ..  1 ) 

CALL  SYMBOL (3 ,0i’? . lLr 2,luY,r . , 1 ) 


rALL  SPALr(  FF,XL“NI,V3,  1 ) 
CALL  £CALE(E;,YLrM,V7, 1) 


. i : r ,a  ' , i i 

CALL  axi«(3,,0.,lH3’  — (Cl«l  .)  ,-11,<lEN,3.,ff<V’*1)  ,Ff<V3«-2)> 

call  axisn.  ,o.,im  , i ,vLen,9u« ,t*(V3+i) , e e ( v y + “* ) ) 

r * i i i T . r ire  f?  . » T 4 n n \ 


CALL  L I P ( r P ,EE,VT,1,  0,0) 
CALL  PLOT  (j», 3) 

PALL  SC  A Lr  ( P''  « XL  CN  , / 1 , 1 » 
CALL  SCAwE  (A  A,  yi.EN  ,'»l  , l) 


SO 


c 'll  axis  n .,  a .,  i?w^ •»--  <ci=-  *.  > ,-i  ?»  xLrr  , o . , *i) , vm  vi  *?n 

CALL  AXISO  . ,0  .,  1 ■*  , 1 , VLEN,90.  , A A ( U*  1) , A A ( Vl*2)  ) 

CALL  SYMBOL  (1 . E,  1 . 2",  . •’fl,  11  , d . ,-l> 

CALL  SYW39L  < .5,  .i  > , . «8 ,11,0.  ,-l> 

CALL  SYMBOL  (1  .25,.  7s?  , . r<8 , 9HULSTA3.  - , C . ,() 

CALL  AFOHD ( 1.  l?v, . d’,  .r;375,  . 96«i~,  . 03  , . Cf  , 15) 

CALL  LINE <99, A A, VI ,1, 0»  0) 

CALL  PLOK3.  ,0., -3) 

CALL  SCALECTC,  YL-:S,'/?,1) 

CALL  SCALE(CC,YLCN,V2, 1) 

CALL  A X 1 5?  ( 3 . ,0  . , 12HC  2-  - (Cl»  - 1 • ) , - l ?,  X LEN.  0 . , 30<  V2*  1)  • Q1  ( V2+2)  ) 
CALL  AXIS(0.,0.,H  , 1 , VL  FU,  9 0 . ,C  ! ( V 2 ♦ 1 ) , Z Z < V 2*? ) ) 

CALL  SYM90w  <1.03i,l.M»,.a»,ll,n.,-l) 

CALL  SVM"OL ( .905, .&3o,  . Ofi  ,11 , G.,-1) 

CALL  SXM9.1L  < 1 .37:  , . rr> , . 0 8 , ^UNSTABLE,  0 3) 

CALL  Af?C(0<  1.23,  .32,  1 1 .,  .03,  . ‘’B,  IS) 

CALL  LlNE<nO  ,CC,V2, 1,0,0) 

CALL  PLOT  ( 0 . -X) 

CALL  SCALE(HL,XLrY,V4, 1) 

CALL  SCALtr(Cr,YL:N>\/,«  , 1) 

CALL  AXIJtO  . ,0  .,  11432-  - <C1*3.>  ,-11 , XLEM,3.  , H*  W,*l  > ,HH  (V-«  *2)  > 
CALL  A XI3 ( 3 • » 0 • » H , 1 , YL  EM, 9 C . ( V,  ♦ 1) , GG  ( V*-  + 2)  ) 

CALL  LINE  (MM, GC,,V4, 1,0.0) 

CALL  PLOTE 
RETUF N 
FNP 
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APPENDIX  F 


Explanations  of  Program 

YFRQOT 

This  program  is  set  up  to  solve  the  polynomial  equation 
which  was  derived  in  the  reduced  linearized  model  analysis. 

It  uses  the  IMSL  subroutine  ZPOLR  which  solves  for  the  roots  of 
a polynomial  equation.  The  accuracy  of  ZPOLR  was  checked  by 
setting  P = 0 and  solving  for  the  roots  of  the  polynomial 
using  the  computer  program  TOTAL.  YFROOT  had  an  accuracy  of 
at  least  10  * which  was  considered  sufficient.  ZPOLR  was 
set  up  to  solve  for  the  roots  repetitively  while  incrementing 
the  roll  rate.  The  inputs  are  specified  in  the  program  YFROOT 
in  Appendix  E.  No  programming  problems  were  encountered. 

EIGEN 

This  program  is,  for  the  most  part,  self  explanatory.  It 
is  limited  to  angles  of  attack  of  0,  10,  20,  and  25  degrees 
and  to  two  altitudes:  0 and  35,000  ft.  Different  altitudes 
can  be  used  if  the  correct  values  of  ROW,  PC,  and  ALT  are  fed 
into  the  program.  EIGEN  uses  the  following  IMSL  subroutines: 
LINV1F,  which  inverts  a matrix;  VMULFF,  which  multiplies  two 
matrices  together;  and  E1GRF,  which  solves  for  the  eigenvalues 
and  eigenvectors  of  a matrix.  Several  test  cases  were  run  with 
simple  2x2  matrices  for  which  the  results  were  known.  This 
was  done  to  verify  that  these  subroutines  did  work.  There 
were  no  programming  problems  except  for  making  sure  that  the 
work  areas  for  the  subroutines  were  dimensioned  properly. 
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This  program  is  set  up  to  plot  four  graphs,  each  with  two 
arrays,  on  the  ASD  Calcomp  Plotter.  The  inputs  are  equations 
which  cannot  be  changed.  All  the  plotting  programs  presented 
the  biggest  programming  problems.  Dimensioning  is  most  important; 
the  array  must  allow  for  all  the  data  points  plus  two  extra 
spaces  needed  for  the  Calcomp  Plotter  routines.  The  largest 
problem  arose  when  trying  to  scale  one  axis  for  more  than 
one  curve.  To  get  around  it,  there  must  be  one  array  which 
contains  all  the  curves  needed  on  one  plot.  This  one  array 
is  used  for  scaling.  Then,  before  calling  the  routine  LINE 
to  plot  these  curves,  the  adjusted  minimum  and  adjusted  delta, 
that  was  obtained  from  scaling  the  larger  array,  must  be  set 
as  the  adjusted  minimum  and  adjusted  delta  for  each  of  the 
curves  to  be  plotted. 

It  is  very  important  to  know  where  the  pen  of  the  plotter 
is  at  all  times.  The  "ASD  Computer  Center  Calcomp  Plotter 
Guide"  is  a useful  aid  for  developing  plotting  programs. 

NONLIN 

This  program  is  set  up  to  solve  for  two  identical  decoupled 
nonlinear  equations.  This  was  done  because  the  subroutine 
NS01A  will  solve  for  a minimum  of  two  equations.  The  inputs 
are  stated  in  the  program.  Again,  dimensioning  is  very  important. 
Arrays  which  are  dimensioned  in  the  main  program  should  be 
dimensioned  the  same  in  subroutine  (s)  that  they  are  used  in. 
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Those  arrays  or  variables  which  are  declared  real,  integer, 
complex,  etc.  must  be  declared  the  same  in  the  subrout ine(s) 
that  they  are  used  in.  In  the  subroutine  PLOTT , the  routine 
AROHD  is  an  auxiliary  routine  of  the  ASD  Calcomp  Plot 
routines.  AROHD  is  available  on  a library  called  CCAUX  and 
must  be  attached  in  the  control  cards. 
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successful,  but  certain  problems  arise  when  working  with  a pair  of 
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of  analysis  known  as  Bifurcation  Analysis  and  Catastrophe  Theory 
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